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THEORETICAL PREDICTION OF PRESSURE DISTRIBUTIONS ON NONLIFTING AIRFOILS AT 

HIGH SUBSONIC SPEEDS ^ 

By JOHN R. SPREITER and ALBERTA ALKSNE 


SUMMARY 

Theoretical 'pressure distributions on nonlijting circular-arc 
airjoils in two-dimensional flows with high subsonic free-stream 
velocity are found by determining approximate solutions, through 
an iteration process, oj an integral equation for transonic flow 
proposed by Oswatitsch. The integral eguation stems directly 
from the small-disturbance theory for transonic flow. This 
method of analysis possesses the advantage of remaining in the 
physical, rather than the hodograph, variables and can be applied 
to airfoils having curved surfaces. After discussion of the deriva- 
tion of the integral eguation and guahtative aspects of the solu- 
tion, results of calculations carried out for circular-arc airfoils in 
flows with free-stream, Mach numbers up to unity are described. 
These results indicate most of the principal phenomena observed 
in experimental studies. Jit subcritical Mach numbers, the pres- 
sure distribution is symmetrical about the midchord position and 
the dray is zero. The magnitude of the pressure coeflicient is 
found to increase more rapidly with increasing Alach number 
than the Prandtl-Glauert rule would indicate. When the critical 
Mach number is exceeded, compression shocks occur, the fore- 
and-aft symmetry of the pressure distribution is lost, and the air- 
foil experiences a drag force. As the Mach number is increased 
further, the shock wave becomes of greater intensity and moves 
rearward alcng the chord, thereby producing a rapid increase in 
the magnitude of the pressure drag coeflicient. jit Mach num- 
bers close to unity, the variation of the pressure, local Mach 
number, and drag conforms, within the limitations of transonic 
small perturbation theory, to the known trends associated with 
the Mach number freeze. Some comparisons with experimental 
results are also included. 

The solutions are obtained using an iteration process which 
differs from the classical methods in that the quadratic nature 
of the integral eguation is recognized. If the iteration calcula- 
tions are started using the linear-theory solution, it is shown that 
the retention of the quadratic feature has the interesting effect of 
forbidding shock-free supercritical second-order solutions. In 
order to obtain solutions for supercritical Mach numbers, it is 
necessary to start the iteration calculations with a velocity or 
pressure distribution which contains a compression shock. When 
this is done, it is found that the iteration procedure converges to 
a definite result. 


INTRODUCTION 

The theoretical problem of transonic flow about thin 
wings has been disciissecl b}^ numerous authors in recent 
A^ears. Since the basic equations are nonlinear and of mixed 
type, the difficulties are great, and progress has been made 
only througli expenditure of considerable effort. At first, 
only tlie basic equations and the similarity rules were estab- 
lished. (See refs. 1 through 12.) More recently, a small 
number of actual solutions have been determined. At the 
present time, the most complete theoretical results are those 
of Guderley and Yoshihara, Vincenti and Wagoner, Cole, and 
Trilling (refs. 13 through 18) for the How about wedge air- 
foils at both subsonic and supersonic speeds. Tliese were 
all obtained by transforming the equations to hodograph 
variables whereby the differential equation becomes linear 
although still of mixed type. Superposition of solutions is 
then possible, but the boundary conditions generally become 
verv complicated. It is liecause of the latter difficulty that 
all the solutions mentioned above are for wedge sections. 
A further disadvantage of the hodograph method is that 
it is definitely restricted to two-dimensional flows, there 
being no known transformation which linearizes the equa- 
tion for three-dimensional compressible flows. 

If the hodograph transformation is not introduced, there 
are available no direct methods of solution. However, 
various iteration methods have been used to study flows 
with high subsonic free-stream velocities. (See ref. 19 for 
a resume.) Almost all these account for the compressibility 
effects by source distributions throughout the flow field and 
start with either the solution for incompressible flow or for 
linearized compressible flow as the first approximation. A 
second approximation is calculated from the first and so on. 
It was not found possible, however, to iterate starting with 
a t}^pical shock-free subsonic flow solution and obtain a 
typical transonic flow field in which the supersonic region 
ends with a shock. 

Oswatitsch has presented another method in references 
20 and 21 for determining the transonic pressure distribution 
on thin airfoils in flows with subsonic free-stream velocity. 
The analysis is carried out in the physical rather than the 


1 Supersedes NACA TN 3090 by John R. Spreiter and Alberta Alksne, March 1954. 
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hoclograpb variables, and leads to a nonlinear integral equa- 
tion in which the unknown velocity appears outside as well 
as inside tlu‘ integral. Oswatitsch finds approximate solu- 
tions not by iteration, but by introducing various functions 
containing undetermined parameters into the integral equa- 
tion and by determining the parameters so that the integral 
equation is satisfied at a small number of points on the air- 
foil. The method is applied to determine the jiressure dis- 
tribution on circular-arc and NACA symmetrical four-digit 
airfoils. The results show certain characteristics of transonic 
fiow such n s the appearance of shock waves and their rearward 
movement across the chord with increasing Mach number. 
It is disconcerting, however, that the method fails to give 
proper res ilts at high subsonic Mach numbers (greater than 
about 0.88 for G-percent-thick circular-arc sections), pro- 
vides a multiplicity of solutions at supercritical Mach num- 
bers, and permits the integral equation to be satisfied at 
only a very limited number of points. 

The present work is based on the integral equation of 
Oswatitsch but an iteration process is used to obtain approx- 
imate solutions. This procedure permits the integral equa- 
tion to be satisfied at a much larger number of points than 
in the original method of Oswatitsch, gives approximate 
solutions 8.t all Mach numbers up to unity, and appears to 
avoid any multiplicity of solutions. The method is applied 
to determine the theoretical ])ressure distribution on sym- 
metrical circular-arc airfoils at zero angle of attack. Excejit 
for phenoriena that are primarily of viscous origin, such as 
boundary- ayer sejiaration, etc., these results exhibit most 
of the experimentally observed features of transonic flows. 

Attention is also called to reference 22 by Gullstrand, an 
associate of (Oswatitsch, in which transonic flows about thin 
air foils are investigated by still another extension of Oswa- 
titsch’s integral-equation theory. Gullstrand sought to 
determine approximate solutions by iteration, although his 
procedure difl'ers considerably from that described herein. 
His method, however, succeeded in determining solutions 
only when tlie Mach number was less than about 0.90 for 
6-])ercent-thick sections. Kesults were given for the pressure 
distribution at sonic speed in a second paper by Gullstrand 
(ref. 23), but they were obtained by introducing a new and 
more comidi(*ated integral equation than that of Oswatitsch 
used herein. In contrast to the present analysis in which 
the entire solution is obtained from the integral equation, 
Gullstrand uses the integral equation to determine only the 
solution for the forward part of the airfoil and uses the 
method of characteristics to complete the solution for the 
rear of the airfoil. Further work of Gullstrand is presented 
in references 24 and 25. 

A list of symbols is contained in Appendix A. 

BASIC EQUATIONS 

The has c equations necessary for the discussion of inviscid 
transonic flow consist of a set of partial differential equations 
relating the velocity components and their gradients at every 
point, together with the auxiliary relation giving the velocity 
jump through a shock wave For thin airfoils inclined at 
zero or small angles of attack, the differential equations can 
be simplifi(‘d by assuming that the shock waves are sufficient- 


ly weak that the flow is irrotational and isentropic, thereby 
permitting the introduction of a velocity potential The 
quasi-linear partial differential equation satisfied by can 
be expressed in the form 

(a^-4>/) 4>..+ (a2-V) 

where the subscript notation is used to indicate differentiation 
and a is the local speed of sound given by the relation 

( 2 ) 

In this latter equation Uo and are, respectively, the free- 
stream velocity and the speed of sound in the free stream, 
and 7 is the ratio of specific heats (for air, 7=1.4). 

It is convenient to introduce the perturbation velocity 
potential </?, where 

iP= — UoX + ^ (3) 

If it is assumed that all perturbation velocities and pertur- 
bation velocity gradients (represented by first and second 
derivatives, respectively, of cp) are small and that only the 
first-order terms in small quantities need be retained, equa- 
tion (1) simplifies to the well-known Prandtl-Glauert equa- 
tion of linear theory 

(1 = ( 4 ) 

where the free-stream velocity is directed along the positive 
X axis as shown in figure 1 and where Mo= Uo/do is the Mach 
number of the free stream. It is well known that equation 
(4) leads to useful results in the study of subsonic and super- 
sonic flow about thin wings and slender bodies but that it is 
incapable, in general, of treating transonic floes. The failure 
of linear theory in the transonic range is evidenced by the 
calculated values of growing to such magnitude that they 
can no longer be j-egarded as small quantities when compared 
with Uo. 



Figure 1. — View of wing and coordinate systerri. 


THEORETICAL PREDICTION OF PRESSURE DISTRIBUTIONS 


ON NONLIFTING AIRFOILS AT HIGH SUBSONIC SPEEDS 3 


vSecond-order tlieory for thin wings would involve solution 
of the equation 

|~<y j ^ ^ 

(1— = I Vz<Pxx+^f^ + 

(5) 


where u, v, and w are the perturbation velocit}- coin])onents 
parallel to the x, y, and 2 axes. This equation corresponds 
to the shock-polar curve for shock waves of small strength 
inclined at any angle between that of normal shock waves 
and that of the Mach lines. On either side of the shock 
wave, the perturbation velocity components are related to 
the perturbation velocity potential in the usual manner 


(See ref. 26, p. 140.) Actually, we are interested in retaining 
higher-order terms only to the extent that is necessary to 
allow the study of transonic flow. Examination of the known 
properties of transonic flow fields indicates that the first 
term on the right can often become of importance and should 
be retained. It is assumed in the small-disturbance theory 
of transonic flow (refs. 1 through 25), however, that the re- 
mainder of the terms on the right can be safely disregarded. 
The simplified equation is 

(1 — ^xx "h " 1 “ ^ 2 Z — ^ ~j r ^x^xx ~ kiPx^xx (6) 

U'o 

where 

(7) 

As a result of minor differences in the perturbation analysis, 
recent papers have used at least three other expressions for 
k. This point will be discussed further at the conclusion of 
the present section. 

Equation (6) is valid only in regions where the necessary 
derivatives exist and are continuous. Since these conditions 
do not hold where shock waves occur, an additional equation 
is needed for the transition through the shock. The neces- 
sary equation is provided by the classical relation for the 
shock polar (e. g., ref. 26, p. 108). 

— ( 8 ) 

7+1 

where i/, r, and w refer to Cartesian velocity components 
with u being parallel to the flow direction ahead of the shock, 
the subscripts a and b refer to conditions ahead of and behind 
the shock, and a* is the critical sound velocity, which can be 
expressed in terms of Uo and Mo as follows: 

jy—] ^ 2 xgN 

C-\7+1+A«7+1) 

The appropriate sim])lified equation is obtained from equa- 
tion (8) by resolving tlie velocities into components parallel 
to the axes of the coordinate system and carrying out a small- 
perturbation analysis analogous to that performed in the 
derivation of equation (6). In this way, the following re- 
lation is found between the perturbation velocity components 
on the two sides of the shock wave: 





( 11 ) 


In addition to satisfying the differential equation and the 
shock-wave equation, the perturbation potential must pro- 
vide flows compatible with the following physical require- 
ments: (a) the perturbation velocities must vanish far 
ahead of the wing and (b) the flow must be tangential to the 
wing surface. Therefore, the following boundary conditions 
are to be specified for the perturbation potential: 


at tlie wing surface W 

Uo\^z)w 


( 12 ) 


(13) 


where bZjbx is tlie local slope of the wing surface in the x 
direction. Furthermore, it is consistent with the assumption 
of small disturbances to satisfy the second boundary condi- 
tion on the two sides of the xy plane rather than on the actual 
wing surface. Equation (13) thus becomes 



where the shape of the wing profile is given by 

c ^ \C s/ 

In addition, it is presumed necessary to prescribe that the 
direct influence of a disturbance in the supersonic region 
proceeds only in the downstream direction and that the Kutta 
condition applies whenever the flow velocity at the trailing 
edge is subsonic. 

Upon solving the above boundary-value problem for the 
potential, one may determine the pressure coefficient for 
planar systems by means of the formula 


r = v—Vo_ 

V p — 

Po TT 2 
2^0 


2 d<p 
Uo c)x 


(16) 


COMPARISON WITH OTHER STATEMENTS OF THE 
TRANSONIC-FLOW EQUATIONS 


(1 -Mo^ {Ua-Uof^ {l'a-V^Y+ {Wa-Wof 

=d// {Uo-Uof 

( 10 ) 


As a result of minor variations in the perturbation analysis, 
recent papers have used at least four different relations for k, 
the coefficient of the nonlinear term in the simplified equation 
for the perturbation velocity potential. As indicated in the 
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preceding- paragrajihs, straightforward development of 
second-order tlieorv loads to the relation 


k^M- 


2 T+I 

U„ 


is sometimes simplified (e. g., refs. 8, 10, and 12) to 


(7) 


(17) 


by arguing that 3/^ can be set equal to unity in this term 
without any loss in accuracy since the right-hand side of 
equation (6) is merel}^ an approximation to allow the treat- 
ment of transonic flows and rapidly diminishes in magnitude 
as Mo de]iarts from unity. In some treatments (e. g., refs. 
7 and 23), equation (1) is divided by a? and the (piotient 
1/a^ in each term is expanded in a binomial series. When 
this is done, the coefficient k of the term involving is 


[: 


2+(7-1)M/ 

Vo 


] 


(18) 


Still auotlier expression for k is used by Oswatitsch in the 
papers that form the princi])al references for the present work. 
Two (leriA'alions are given, one based on mass-flow consider- 
ations (ref. 21) and the other (ref. 20) on simplifying equation 
(1) under the assumption of nearly ])arallel flow to 


(1 — ^P)ipxx'\~ ^vy-\~ 


(19) 


and substituting the following series for the variable coeffi- 
cient (1 — 3/^): 


1 -3/2=1- 


^3/,2 
a"* 


^ f2 t A ® I 


( 20 ) 


where 3/ is the local Mach number and a* is the critical 
sound velocity as defined in equation (9). Comparison of 
equations (19) and (20) with equation (0) shows that the 
coefficient k in this approximation is 


1-34 2 


(21) 


through 6, 9, and 13 through 18.) It is shown in reference 
12, however, that tlie pressure, force, and moment results 
ol)tained using tliese equations are identical with those 
obtained using the present equations if k is selected as given 
in equation (17). These results, however, can be easily 
converted to those that would be obtained using any of the 
other expressions for k by simply re])lacing {y+l)/Uo by k 
wherever it occurs. 

In order to facilitate comparison with previous results and 
to achieve an economy of notation, the present anal^^sis is 
carried as far as possible without specifying a particular 
relation for k. That is, the equations of the analysis and the 
reduced parameters with which the results are expressed are 
written containing k which may be equated to any of tlie 
four stated exj)ressions. However, the actual values of the 
])ressure coellicient and Afacli number for an airfoil of specific 
thickness ratio depend on which relation is selected for k. 
Idle present calculations have been made using the expression 
for k given in equation (7). The principal reason for this 
choice is that it appears to provide a set of equations, or a 
mathematical model, which approximates certain essential 
features of transonic flow with superior accuracy. I^eforc 
proceeding with discussion of this point, it should be noted 
that the four alternative expressions for k are identical for 
Mo=lj and all but that given by equation (17) are zero for 

34 = 0 . 

A significant case where the four relations lead to different 
results is the prediction of the variation with free-stream 
Mach number of the critical pressure coefficient defined 
as the value of the pressure coefficient Cp at a point where 
the local Mach number is unity. It is important that a 
reasonably good approximation be maintained for the varia- 
tion of Cp^^ with 34 because shock waves make their first 
appearance and the airfoil first experinces a pressure drag 
when Cp becomes more negative than Cp^^ somewhere on the 
airfoil surface. In the present approximation, Cp^^ corre- 
sponds to that value of Cp and, hence, of <^-, at which equa- 
tion (6) cfianges locally from elliptic to hyperiiolic type. 
This condition is recognized by the vanishing of the coefficient 
of ipxxy thus 

l-3//-A’(^x)c.=0 


A similar situation arises in the derivation of the simplified 
equation for the shock ])olar. Here again the precise form 
of the expi-ession for the coefficient k of equation (10) depends 
on the details of the perturbation analysis. The most 
important point from the present point of view is that the 
same expression for k is used in both the equation for the 
potential and that for the shock polar, namely, equations 
(0) and (19). AVhile this point has not always been explicitly 
stated, it is actually a necessary condition for tlie existence 
of the well-known transonic similarity rules. 

The foregoing discussion has been based on equations 
obtained by assuming that the local velocities are only 
slightly diff(U*ent from the free-stream velocity. On tlie 
other hand, many of the recent papers on transonic flow 
about wings and bodies have been based on equations ob- 
tained b}^ assuming that the local velocities are only slightly 
different from the critical sound velocity a*. (See refs. 1 


or, in view of equation (16) 



(<^x)cr = — (1—342) 


( 22 ) 


The exact relation for isentropic flow is (e. g., ref. 27, ]). 28) 


C 

1) 



7—1 

7+1 



(23) 


The variation of Cp^^ with 34 has been computed using the 
exact relation and each of the four approximate relations. 
The results are ])resented graphically in figure 2. It may 
be seen that a reasonably good approximation for Cp^^ is 
obtained over a wide Mach number range when k is taken 
as 3/o2(7+l)/f7^ or (1 — 34^)/(«*— f4), and that a some- 
what greater error is incurred when k is equated to 3//[2 + 
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Figure 2. — Variation, approximate and exact, of critical pressure 
coefficient with free-stream Mach number. 


{y—\)Mo^]IUo- On tlio otlier liand, a very poor approxi- 
mation results if k is equated to (y+O/t^o- 

Similar comparisons can l)e made for local Mach numliers 
M otlier than unit}' by noting that the coeflicieiit (1 — 
]\l^—k(px) of (pxx in equation (6) corresponds, in the present 
approximation, to 1— dP, thus, 

hTT 

C, (24) 


The corresponding exact relation for isentropic flow is 




r y—\ > 

7 ” 

7-1 

-1 + 

z 




_ 

L 2 J 

_ 


(25) 


The results so obtained are generally similar to those indi- 
cated in Figure 2 although the relative accuracy of the 
better approximations changes somewhat with the situation. 
All the approximations are exact, of course, when 
On the other hand, none of the approximations are exact, 
except for isolated cases, when Cp is different from zero, 
even though all of the ap|)roximations agree among them- 
selves when the free-stream Mach number is unity. In 
order to provide some information regarding the errors that 
are likely to be incurred when Cp is not very small, figure 3 
has been prepared illustrating the variation of local Mach 



Cp 

Figure 3. — Variation, approximate and exact, of local Mach number 
with pressure coefficient, iV/o=l. 


number with pressure coefficient for a free-stream Mach 
number of unity. 

A second case where the exact and approximate relations 
can be compared is furnished by considering the velocity 
jump through a shock wave. If the flow ahead of the shock 
wave is uniform and pai*allel to the x axis, the results may 
conveniently be represented by the shock-polar diagram in 
which is plotted as a function of The exact 

relation is furnished by equation (8). The corresponding 
approximate relations are determined from equation (10) by 
setting Ua, Va, hiid \Oa to zero, whereby 

w,] (26) 


Once the variation of with Ujj is determined for a 

given Mo, the corresponding variation of with u* 

may be readily determined since, for this case, 

u^=Uo+Ut, Vb=i% Wb=Wb (27) 


The variation of with Ui, for Mq=1-2 has been 

computed using both the exact and approximate relations, 
and the results are presented graphically in conventional 
shock-polar form in figure 4. 

It is evident from this comparison that the best approxima- 
tion to the shock-polar curve is that obtained by equating 
k to Mo^{y+l)/Uo- Since all shock waves are assumed to 
be normal to the flow direction in the course of the present 
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Figure 4. — Approximate and exact shock polar diagram, A/o=1.2. 
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analysis, a notable point is that this expression for k leads 
to the exap.t relation for the velocity jump through a normal 
shock wave. 

In problems such as we are considering here, the final test 
is providerl by comparison with experimental results. Al- 
though both experimental and theoretical results for the 
transonic speed range are limited, complete information 
does exist at the present time for the drag of a single-wedge 
section followed by a straiglit section extending far down- 
stream. (See refs. 13, 14, 16, 28, and 29.) Tlie theoretical 
results we]*e determined originally using ec|uations obtained 
by assuming that the local velocities are only slightly 
different from the critical sound velocity a* and are therefore 
identical \\ ith those that would be obtained using the present 
eciuations, provided k is equated to ( 7 +l)f^o. Figure 5 
shows the theoretical and experimental results plott(‘d in 



Figure 5. — 'Pheoretical and experimental drag results for a single wedge 

7 + 1 

section, fc=— * 

the same manner as in the original papers. The small 
vertical lines on the experimental data points represent the 
uncertainty <^f the values. This figure indicates that the 
theoretical and experimental results are only in general 
qualitative agreement when k={y-\-\)Uo. 

Tlie same results are replotted in figure 6 with k equated to 
il//( 7 +l)/r/„ rather than ( 7 +l)/f^o. It can be seen that 
the theoretical and experimental results are now in nearly 
perfect agreement. Comparison of figures 5 and 6 provides 
striking evidence supporting the contention that k should be 
equated to d/o^( 7 + 1 )/f/o rather than ( 7 +l)/f/o. 

DERIVATION OF INTEGRAL EQUATIONS FOR TRANSONIC 

FLOW 

In order to make the present work more self-contained, a 
derivation of the integral equations for transonic flows 
having sulisonic free-stream velocities will be presented even 
though this has been done previously by both Oswatitsch 
and Gullstj-and (refs. 20, 21, and 22). The present deriva- 
tion, in common with that of Gullstrand, proceeds through 
the application of Green’s theorem in a manner closely 
analogous 1o that emplo 3 ^ed in linearized wing theory (e. g.. 



Mq^-\ 

[(y + l)/W/(//cy^ 


FTgure 0. — Theoretical and experimental drag results for a single wedge 

; A/oHT + I) 
section, k= ^ 


ref. 30), except that proper cognizance must be taken of the 
shock discontinuities and of the additional nonlinear term in 
the differential equation for the perturbation velocity po- 
tential. For the sake of completeness and to illustrate the 
simplifications introduced by making additional restrictions 
and assumptions, the derviation will be carried througli for 
lifting wings of finite thickness and span, even though the 
applications contained herein will be confined to two- 
dimensional flow about symmetrical airfoil sections at zero 
angle of attack. 

The differential equation fundamental to the following 
discussion is equation (6). 






( 6 ) 


Siiice the principal object of the following analv'sis is to 
determine the pressure which, according to equation (16), is 
linearlv proportional to the pert uihat ion velocity component 
it is (convenient to work with an (equivalent equation for u 
obtained b\' differentiating equation (6) with respect to x; it is 




(f) 


(28) 


It is advantageous to noiinalize tlie e(|uations by letting 
x=x, y=^y, z=0Z: 


_ A- 


- d(f k 


dx 


_ ^ip k 




where 


In this wa\', e(iuatioiis (6) and (28) ]*educe to the following: 

(30) 
5“ 


^ ^ 2-^^ dV 


d'u , , d'u 

dy 


d'-u . d^u . d^u 5“ /u \ 


(31) 
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Before proceeding, it should be noted that the introduction 
of the reduced perturbation velocity component u permits 
the ready recognition of regions of subsonic and supersonic 
velocities and emphasizes the points at which sonic velocity 
occurs. This relationship becomes immediately apparent 
upon substituting the definition of u into equation (24). 
Thus, 


1 


-Mo 


, u=l — u 


(32) 


from which it is clear, for flows having subsonic free-stream 
Mach numbers (Mo<^l), that u<Cl when the local velocity is 
subsonic, u=l when it is sonic, and u^l when it is super- 
sonic. 

As noted above, the derivation proceeds from Greenes 
theorem which relates a volume integral over a region B to a 
surface integral over the surface 2 enclosing B. If o- and 
are any two functions which, together with their first and 
second derivatives, are finite and single-valued throughout 
B, Greeids theorem states 



dn dn / 



fj'l' dli 


(3:^) 



Ficure 7 . — Region of integration. 


where the directional derivatives on the left side are taken 
along the normal n, drawn inward, to the surface 2. It is 
convenient to let ^=u and to choose a as the fundamental 
solution l/rs of the equation W = 0 


1 1 

(jj-vr+ 


(34) 


whereby equation (33) becomes 



a m-i 

dn dn Vra/j 


dX=r= — 



/.fj 


(35) 


The variables of integration in the equation are f, v, f while 
S, i are the coordinates of a point P. It must be observed 
that l/ra is singular at ra=0 and u is discontinuous at the 
shock wave. The point P and the shock wave must, there- 
fore, be excluded from the region B. llie exclusion of P 
from the region B is accomplished b}^ enclosing it within a 
small sphere. The shock wave is excluded by altering the 
boundary of the region so that it goes around the shock wave. 
In this way, equation (35) may be applied to the region Bu 
bounded by the xy plane and a hemispherical dome of infinite 
radius bung above this plane, exclusive of the subregions 
surrounding P and the shock wave (see fig. 7). Since, 
furthermore, the values of u may be assumed to diminish 
sufhciently rapidly with distance that the contributions of 


the integrals over the hemisphere vanish, the following 
result is obtained: 



where the subscript u denotes conditions on the upper side 
of the xy plane, the subscripts a and h denote values immedi- 
ately ahead of and behind the shock wave, and S is the sur- 
face of the shock wave. The volume integral is defined as 


follows when P is ahead of S. (^For sake of brevity, is 


written in place of 


i 

rzdf 


(I)') 


Hu 

J ^co ■ r*Xp-tx . C^s -^1 - A 

\ ( ^d^-{- I |_ ^d^j (37) 

— 00 •/ 0 d — 00 


It is clear that the corresponding definition of the volume 
integral when P is behind S can be obtained by rearranging 
the limits of integration. If P is kept fixed in the upper 
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half space and tlie region Ri bounded l)y the xy ])lane and a 
liemis])heric.il dome of infinite radius lying below this plane 
is considered, it follows in a similar manner that 


45rJJ l_r3 

— 00 





(38) 


where the subscript I denotes conditions on the lower side 
of the xy plain' and the volume integral is defined as follows: 


lU ^ lU 


(lit ions; for example, (a) given the loading on a lifting sur- 
face, find the camber distribution; (b) given the shape of a 
symmetrical nonlifting airfoil, find the pressure distribution. 

Ecpiation (41) can be regarded as the final integral equa- 
tion for t/, but it is advantageous for the forthcoming anal- 
ysis to ])erform two more operations. They are to integrate 
the volume intc'gral twice by parts with respect to i, taking 
proper cognizance of the definitions given in equations (37) 
and (39), and to decompose the surface integral over the 
shock wave into components parallel to the axes of the 
coordinate system. In this wa}^, the following equation is 
obtained: 


'“c JJ[k (r.)] 


J +» ro . f rxs-<2 - r+“ -'i 

<^\ 1 Hi y 

-CO •/—CO CD «/ « 2 J 


Introducing: (lie notation 


. bu dUu bill 
bu=u,—vi, A — 

bi bt bt: 


(39) 


(40) 


and addinjr cainations (30) and (38), we have 


(41) 


where the surface integral over S represents the sum of the 
corresponding integrals over Su and Si and the volume 
integral over Ji represents the sum of the integrals over Ji\ 
and Hi. The integrand of the first integral of equation (41) 
is z(To over all of the xy ])lane except the surface W of the 
wing sufficicuiily extended to include the effect of the edge 
singularities and is, in certain cases, exactly ecpial to the 
value of u given by the linearized theory of subsonic flow 
about thin wings (e. g., ref. 30). 



It can be S(‘en that the first integral of equation (41) may 
be equated o when the problem is one in w^hich 

Au = AUl 

and 

Mdu/dJ) = A(du/df)L 

'’riiis condition exists in those problems where Au and 
A(dii/df) are jirescribed at the outset by the boundary con- 



Although the integi’ation by j)arts of the triple integral 
])erformed in going from equation (41) to (43) may seem 
somewhat arliitrary, the resulting equation is superior from 
the point of view of obtaining approximate solutions. For 
('xample, the triple integral of ec] nation (41) shows a veiy 
strong influence of the velocities in the region immediateh^ 
surrounding P since the}- are multiplied by l/vs. This in- 
fluence is largely nullified in the triple integral of equation 
(43) because part of the region has a negative influence and 
part has a positive influence. The predominant influence 
in the latter case is furnished by the term u~l2 standing out- 
side the integral. The contribution of distant regions is 
also diminished in importance in the trijfle integral of equa- 
tion (43) since their influence varies inverseh' with the third 
power of the distance, rather than the first power as in 
ectuation (41). The advantages of using the formulation 
])rovided by equation (43) will become more evident on 
examining A])pendix B and the section entitled “Sub- 
critical flows. 

A further advantage is that the value of the triple integral 
of eciuation (43) is continuous through a shock wave rather 
than discontinuous as is the case with equation (41). A 
point of great importance in the approximate solution de- 
scribed herein arises from the fact that the integration by 
parts provides extra terms (those containing IP 12) in the 
integrals along the shock surface S which combine with those 
already present in such a way that the contribution of these 
integrals becomes very small when the shock waves approach 
normal waves, as is iisualW the case at high subsonic speeds. 
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In addition to satisfying the integral equation for u given wliere 
in equation (43), the velocity components on opposite sides 
of shock waves must be in accord with the simplified relation 
for the shock polar given in equation (10). This equation and 
may be rewritten in normalized form by introducing the 
quantities defined in equation (29), thus 


1 _ 1 


',49) 


(50) 


{Ua-U,f (44) 


Two alternative forms of equation (44) arc the following: 
and 

{Ua—U^ [^(wa— ^ — (tii.— ^)J+(F„— »6)^+(w„— Ws)^=0 (46) 

If the shock wave is a normal wave and the fiow is parallel 
to the X axis (i.e., Va=Vb=Wa=Wb=^, hut Ua9^Ub)i it can be 
seen from equation (45) that the normalized perturbation 
velocity component u jumps from 1+A immediately ahead 
of the shock to 1 — A immediately behind the shock. On 
the other hand, equation (46) shows that the quantity 
is equal on the two sides of the shock. This is consistent 
with the fact that the latter quantity corresponds, in the 
transonic approximation, to the mass flow, which is continu- 
ous through a normal shock. 

The solution, by the present method, of the general prob- 
lem of transonic flow about thin wings requires the solution 
of equation (43) while taking ])roper account of the shock 
relations given in equation (44). This represents a formida- 
ble task well beyond the reach of the present analysis. 
Simplification can be achieved in two ways: by restricting 
attention to a less general class of problems and by introduc- 
ing additional simplifying approximations. The first way is, 
of course, much to be preferred. Accordingly, in most of 
the following analysis, attention will be confined to two- 
dimensional flows. The necessary equations can be obtained 
from equations (42) through (44) above by integrating in the 
direction from to noting that ^=0 and 

that u and w are independent of ry. They are as follows: 


“s If I? I)''* 
i JJI 5- s { [“' I 

R 

V 2/5J rjo LV bi r2/a 

In i I, In i) 1 f \ rff (48) 

^2/6j Lcos (?^,?)J* j 


As remarked following equation (42), considerable sim- 
plification results in both two- and three-dimensional prob- 
lems if attention is confined to the determination of the pres- 
sure distribution on symmetrical noidifting wings of specified 
geometry. This restriction permits the introduction of the 
relations 

AiZ=Ai7z,=0 





(51) 


into the integral over W in either equation (43) or (48). 
This integral is then equal to the linear-theory solution Uj^ 
given in equation (42) or (47) and can be determined com- 
pletely at the outset of the analysis. 


SIMPLIFICATION AND APPROXIMATE SOLUTION OF 

INTEGRAL EQUATIONS FOR TRANSONIC FLOW ABOUT 

NONLIFTING WINGS 

QUALITATIVE DISCUSSION OF INTEGRAL EQUATIONS 

The integral equations and the auxiliary relations develop- 
ed in the jireceding section provide a means for the discus- 
sion of the aerodynamics of symmetrical nonlifting thin 
wings of specified geometry in flows with free-stream Mach 
numbers up to unity. At the present stage of development, 
however, it is necessary to introduce some further approxi- 
mations before solutions can be obtained. One of the more 
accurate of tliese involves two statements about the nature 
of the shock waves. They are: (a) all shock waves are 
assumed to lie in a plane perpendicular to the x axis, and 
(b) the shock waves are assumed to be normal shock waves 
(i. e., normal to the local flow direction). These two state- 
ments are slightly contradictory in themselves but might 
be expected to approach the true conditions quite closely 
for flows about smoothly curved thin airfoils. The first 
statement corresponds to setting cos(/i, ^) and cos (?i, J) to 
zero, thereby eliminating part of the integrals over S of 
equations (43) and (48). The second permits an advan- 
tageous introduction of equations (30) and (46) to eliminate 
the remainder of the integrals over S. The above assump- 
tions correspond to setting v and w to zero both before and 
after the shock wave and lead to the following relations: 


Equations ('43) and (48) thereby simplify to 

“==‘+¥-4;XfjT (I) 

R 

for three-dimensional flows, and 


(52) 


(53) 


( 54 ) 
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for two-dimensional flows. Those equations correspond to 
those used by Oswatitsch and Giillstrand (refs. 20, 21, and 
22) although a niiinber of further assumptions were necessary 
before approximate solutions could he obtained for the 
velocity distribution on an airfoil surface. The present 
analysis also requires many of the same or similar assump- 
tions, but there are a number of general points which should 
be discussed before any further approximations are intro- 
duced. Th(‘se points were not observed in the previous 
papers and tlie results suffer by lacking uniqueness in some 
cases and failing to converge in other cases. 

Before proceeding further, it should be observed that the 
solutions of equations (53) and (54) must approach those of 
linear theory when the free-stream Alach number is much 
less than unity, since u<^\ and the terms involving the 
square of u become negligible with respect to those linear in 
Uj thereb}^ leaving only 




(55) 


In the interest of brevity, the integral equations for u 
given in equations (53) and (54) may be rewritten as follows: 

Tr 

/ 

(56) 

'2 

where 

i)4<wl 

(57) 

for three-dimensional flows, and 




(i)'«] 

(58) 


for two-dimensional flows. Although / is a function of u 
and is thei-efore unknown, it is informative to rewrite 
equation (50) by solving for u in terms of / and ul, thus 


u = \±^I-{2u^-\) = \±^I-L (59) 


Several points are to be observed at once with regard to 
equation (59). First of all, the discriminant must always 
be positive in order to obtain real values for u, thus 

1>L (60) 


Furthermore, the choice of the plus or minus sign determines 
whether the local velocities are subsonic or supersonic. A 
change in sign at a point where the radical is zero corresponds 
to a smooth transition through sonic velocity. A change in 
sign at a point where the radical is not zero corresponds to a 
discontinuous jump in velocity. As pointed out following 
equation (46), such discontinuities correspond to normal 
shock waves and are permissible when they proceed from 
supersonic to subsonic velocities (or from plus to minus sign 
in equation (59)) when progressing in the ffow direction. 
Discontinuities in the reverse direction are inadmissible since 
they correspond to expansion shocks, a phenomenon which 
violates the second law of thermodynamics. 


The values of and hence L, can be calculated for any 
given wing and are generally characterized by certain regions 
in which is positive and other regions in which it is nega- 
tive. The absolute values increase continuously with in- 
creasing Mach number and the maximum positive values 
may considerably exceed unity as sonic velocity is approached 
ill the free stream. Not very much can be stated at this 
point al)out the values of /, except that they dejiend on the 
distriliution as well as magnitude of u and that the above 
inequality must be satisfied. The relation between the two 
curves is of utmost importance, however, and will he dis- 
cussed qualitatively in the following paragraphs. 

In order to remove unnecessary complications and to facili- 
tate the discussion, the following remarks will lie confined to 
the relations in two-dimensional flow lietween the functions 
I, L, u, and evaluated at the airfoil surface. In this way, 
each of the four functions reduces to a function of a single 
variable x and can be illustrated simply by curves rather 
than surfaces or hypersufaces. The subscript IF is appended 
to u and Ul to denote that the values are those at the wing 
surface. In that which follows, the curves will be shown on 
two separate jilots, one containing the L and I curves rep- 
resenting the components involved in the solution of equa- 
tion (59) and the other containing the u^r and curves 
representing, respectively, the velocity distributions given 
by transonic theory and by linear theory. In order to make 
the discussion more definite, the curves will be drawn quali- 
tatively as they would appear for a circular-arc airfoil having 
its maximum thickness at the midchord position. A quanti- 
tative discussion of these characteristics will be taken up for 
the same airfoil in a later section. 

The linear-theory solution for subsonic ffow about cir- 
cular-arc airfoils can lie easily derived through application 
of the expression given in equation (47). It is found that 
the values of are symmetrical about the midchord posi- 
tion at all free-stream Mach numbers less than unity. It 
follows directly that the L curves possess the same symmetry. 

For pure subsonic flows about such an airfoil, it is well 
known that the more exact treatments, such as the Karman- 
Tsien or the classical iteration methods reviewed in reference 
19, indicate that the nonlinear-theory solutions and, hence, 
the uw and / curves, are also symmetrical about the mid- 
chord position. Since sonic velocity is not attained at any 
point, the L and 1 curves never touch. Sketches of the 
curves for this condition are shown in figure 8. The second 
part of this figure illustrates the fact that uw possesses larger 
values over the middle of the airfoil chord than does UL^y, 




ii, and ul in the subcritical range. 
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This relation has not been dednced from the foregoing pre- 
liminary considerations but is known from inaiu^ sources in 
classical subsonic theory. The same result has been found 
once again in the present work by canying out an approx- 
imate solution of equation (59). The details of these calcu- 
lations will be described in a later section. 

The curves shown in figure 8 are typical of those for all 
Mach numbers less than the critical Mach number 
defined as the lowest free-stream Mach number at which 
sonic velocity (Jiw=l) occurs somewhere on the wing surface. 
The aforementioned symmetry properties are preserved for 
all ^^ach numbers less than the critical, but the amplitudes 
of all four curves increase with increasing Mach number. 
For the symmetrical circular-arc airfoils considered here, the 
maximum values of and uw and, hence, L and / occur at 
the 50-percent-chord position for all subcritical Mach num- 
bers. Thus, in addition to the requirement that I>L at 
every point, it is necessary when the Mach number is equal 
to the critical value, that I=L at the 50-percent-chord posi- 
tion. A rather interesting additional requirement that fol- 
lows from the quadratic nature of equation (56) together 
with the assumption that u reaches a smooth maximum at 
the 50-percent-chord station is that the / and L curves have 
not only the same first derivatives at this station but also 
identical second derivatives. A typical set of curves for this 
Mach number is shown in figure 9. It is to be remarked 
that the minus sign is to be used in equation (59) for all 
free-stream Mach numbers equal to or less than M^r- 




Figure 9. — Typical curves of /, L, u, and ul at the critical Mach 

number. 

It is interesting to contemplate the various possibilities 
that may occur when the Mach number is increased beyond 
the critical. Accordingly, let us first consider the implica- 
tions of assuming that the velocity distribution remains 
symmetrical and continuous across the chord and simply 
increases in magnitude with increasing Mach number. 
(This is, in fact, the only type of mixed flow field that the 
classical iteration methods have indicated, but many 
doubts have been expressed from time to time regarding 
the convergence of the methods when Mo^Mcr-) ith such 
an assumed symmetry of the flow, the four curves will 
appear qualitativeh" as shown in figure 10. As may be 
seen, the curves are all rather similar to those previously 
discussed for lower Mach numbers. The outstanding 
difterence is the relation between I and L. When Mo=Mcr, 
the two curves are tangent at the 50-percent-chord station, 
and the radius of curvature of the / curve at the same sta- 




Figure 10. — Typical curves of /, L, w, and ul corresponding to a 
shock-free supercritical flow, symmeterical solution. 

tion is equal to that of the L curve. When 
tangency occurs at two points, equally spaced before and 
after the 50-percent-chord station, and the sign in equation 
(59) is to be changed to plus over the portion of the chord 
lying between the two points of tangency. In order for 
Uw to attain its maximum value at midchord, as shown in 
figure 10, the radius of curvature of the I curve must be 
less than that of the L curve at the 50-percent-chord station. 
At the tangent or sonic points, it follows from equation 
(59) and the assumption of smooth acceleration or decelera- 
tion through sonic velocity that the second derivative of 
the / curve is greater than that of the L curve. Similar 
considerations apply for airfoils that are not symmetrical 
about the midchord station. The occurrence of shock-free 
flow would again require that the / curve be tangent to the 
L curve at two points along the chord. 

Before leaving this subject, it is interesting to inquire 
what the result would be of a slight alteration of the airfoil 
shape, assuming that the original shape was sucli that the 
associated flow was of the shock-free mixed type. Consider, 
for sake of simplicity, that the basic airfoil is symmetrical 
about the 50-percent-chord station. The assumed shock- 
free flow is therefore symmetrical about the midchord 
station and appears as shown in figure 11. In order to 
preserve the geometrical symmetry, consider that the air- 
foil shape is changed by the addition to both upper and lower 
surfaces of small bumps located at the 50-percent-chord 
station. Since the flow adjacent to the bumps is supersonic, 
their disturbance pattern is propagated downstream in nar- 
row bands which reflect alternately from the sonic line 
and the airfoil surface as shown in figure 12. Busemann 



\ / 
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Figure 11. — Shock-free transonic flow. 
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Figure 12. — Transonic flow with shock wave. 

has pointed out in reference 31 that as the disturbances 
reach the I’car of the supersonic region where the local 
Mach numbers of tlie original How approach unity, the 
reflections become increasing!}^ concentrated, the disturb- 
ances aiipdify, and finally a shock occurs. It thus appears 
that there* are only a restricted number of shock-free mixed 
flows and that, in general, the supersonic region terminates 
with a shock wave. 

On the other hand, if the disturbances were projiagated 
forward they would culminate in a shock wave at the forward 
sonic point. This shock, however, would be an expansion 
shock in which the velocity jumps from subsonic to super- 
sonic values. As mentioned previously, expansion shocks 
are forbidden by thermodynamic considerations and, hence, 
must be (‘xcluded from the present analysis. This can be 
accomplis.ied by stipulating that the transition from sub- 
sonic to supersonic velocities at the forward sonic j)oint be 
smooth, or that the L and / curves retain one point of 
tangency. In addition to preventing the occurrence of ex- 
pansion shocks, it appears, on the basis of the foregoing, that 
this requirement also effectively rules out all the undesired 
forward propagating disturbances. 

The preceding discussion provides an insight into the 
mechanism for the develojiment of asymmetrical flows about 
symmetrical airfoils and for the occurrence of shock waves. 
In the present work, these considerations are reflected in the 
relation b etween the / and the L curves. Thus, recall that 
if the flow is shock free as shown in figure 11, the / and the Z 
curves an* both symmetrical as shown previously in figure 
10. If th ^ airfoil shape is changed in the manner indicated 
in figure 12, it is evident that both the / and the Z curves 
will become altered. In so doing, the Z curve remains sym- 
metrical about the midchord statioii, but the asymmetrical 
nature of the u distribution causes the I curve to take on 
larger values over the rear of the airfoil than over the front. 
If the flow adjusts itself so that the / curve is tangent to the 
Z curve at a point on the forward half of the airfoil, as it 
must do to avoid the occurrence of forbidden expansion 
shock w ^A es, it is likely that the curves w ill not be tangent 
at a second point along the rear half of the airfoil. In the 
application of equation (59), therefore, the sign changes 
from mini s to plus at the point of tangency but may change 



back to minus at a point wdiere the tw o curves have difterent 
values. As noted previously, such a condition corresponds, 
in the present analysis, to the occurrence of a shock w’ave 
and is associated with a curve of the type shown in 
figure 13. 

Since it has been indicated that the shock-free mixed flow 
is the exception rather than the rule, the curves of figure 13 
may be regarded as typical for the supercritical Mach num- 
ber range. In any event, these curves may lie considered 
as the more general ones since they include the syminetrical 
curves as a special case. These matters will arise again and 
be the topic for further discussion in the next section in which 
an apjiroximate method for the solution of the integral 
equation for transonic flow' w ill be described. 




Figure 13. — Typical curves of /, L, u, and Ul corresponding to super- 
critical flow^ with shock, asymmetrical solution. 

SIMPLIFICATION OF THE INTEGRAL EQUATION 

The remainder of the present discussion w ill be concerned 
wdth two-dimensional flow about nonlifting symmetrical 
airfoils of specified geometry under the assumjition that any 
shock w aves which may be present are normal shocks situated 
perpendicular to the x axis. The analysis w'ill be based there- 
fore on equation (54) w'hich, w'hen w ritten in full, is 


-|- 00 -j- 00 



00 — 00 
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(61) 


Approximate solutions of this equation could conceivably be 
wmrked out numerically by starting with a tw'o-dimeiisional 
grid of suitably selected values for u(^,f) and iterating until 
convergence is obtained. Such calculations w'ould proceed 
])V inserting the assumed values for u into the double integral 
and solving to obtain the next aiiproximation for u{x,z)f 
making use of the tangency condition on the surfaces or 
functions represented by / and Z as discussed in the preceding 
section. If the first approximation for u is taken to be the 
results given by incompressible or by linearized compressible 
flow' theory, as in the Kayleigh-Janzen and other classical 
iteration methods, it seems to be the prevailing belief that 
convergence will be obtained only wdien the free-stream 
Mach number is sufficiently small that the flow' is subsonic at 
every point. It is at this point that Osw'atitsch (refs. 20 and 
21) supplied the important idea that mixed flow' fields con- 
taining shock w'aves can be obtained if the starting u distri- 
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biitioTi contains shock waves. Tims, in place of taking the 
starting solution to be that given by incompressible or 
linearized compressible flow theory, the idea is to start with 
a reasonable guess for the values of u, being sure to include a 
proper discontinuity complying with the shock relations of 
equation (52), and then to ])rocced to the solution. As will 
be seen in the succeeding sections, it is not necessary to be 
highly accurate in the initial guess for u. 

A source of difficulty in the numerical solution of equation 
(61) by an iteration process is the double integral. If it 
could be reduced to a single integral by introducing a suitable 
approximation, the entire problem would be greatly simpli- 
fied. In the present analysis, it is assumed, following 
Oswatitsch, that approximate knowledge of the velocity 
distribution is sufficient for providing a working approxima- 
tion for the double integral. In particular, it is assumed that 
a sufficiently good approximation to the velocities in the 
vicinity of the wing can be expressed in terms of the local 
coordinate i, the ordinates of the airfoil surface i(x), and the 
desired but unknown velocity distribution Uw(x) on the 
airfoil surface. This permits one integration to be performed, 
thereby reducing the double integral of equation (61) to a 
single integral. 

A number of statements regarding the variation of u with 
i over the middle portion of the airfoil can be made immedi- 
ately. For example, u starts from the value at the air- 
foil surface with an intitial rate of change given by the 
irrotationality condition 

(iXKil 


and probably vanishes at great distances as 1/i^. These 
conditions, of course, are not sufficient to determine com- 
pletely the variation of u with but may be used as the 
basis for the development of an approximate relation. 
Oswatitsch (refs. 20 and 21) has already considered this step 
and has suggested the following relation : 


i/(x,i) 


Uw(x^d) 

[^+mY 




where b is a function of x so chosen that the irrotationality 
condition is fulfilled at i=0. Thus, differentiation of equa- 
tion (63) with respect to i and insertion of the definitions of 
equation (29) and the boundary condition of er{uation (14) 
yield the following: 

7 ^Uw U\v 

(di/di)ir (dlc/dir)^- k {6wldx)w 

2u^ . . 

kUo &Zjbx^) (yz/dP) ^ ^ 

where Z represents the reduced ordinates of the airfoil, 
related to the actual ordinates by 

Zm=-^Z(x) (65) 


Attention is called to the fact that the approximate relation 
for u(x,z) given above is not entirely satisfactory. Evidence 


of this is provided by the fact that u is indicated to be zero 
in the region ahead of the leading edge and behind the 
trailing edge where b is infinite and that the discontinuities 
in u at the shock surfaces are consistent with the shock 
relations only at the surface of the airfoil. The errors in the 
pressures on the airfoil surface resulting from the former are 
small due to the attenuating influence of distance, and those 
resulting from the latter have been partially compensated 
for by a readjustment of the approximation at Mach num- 
bers near unity where the shock strength becomes large. 

Substitution of equation (63) into the double integral of 
equation (61) permits integration with respect to f. Thus, 
by performing this integration and setting z = 0, the following 
approximate integral equation is obtained for 

(66) 

The function E is 

E [I |X|( 5 -. 0 ^-+X,- 

(1 - 10 .V+oZ^)ln |V|-L ( 1 +X 2 )( 25 - 71 Z^-X^-X«)J 

( 67 ) 

The nature of E(X) is illustrated graphically in figure 14. 



A 

Figure 14. — Variation of E{X) with X. 
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Although the integration interval is indicated in equation (61 ) 
to extend from J = — oo to^ = +oo^ the contribution of the 
regions ahead of the leading edge and behind the trailing edge 
is zero sinc(‘ h is infinite. The integration need, therefore, be 
carried out only over the chord c. It should be noted that 
the integral in equation (66) corresponds to //2 in equation 
(56). 

Although e((uation (6()) is considerably sini])Ier than 
e({uation (61) owing to the re])lacement of the double integral 
by a single integral, many of the essential difliculties remain 
since the integral ecpiation is still nonlinear and the kernel 
is singular. Since no known methods exist for the solution 
of such etpiations, we can only j)roceed at the j)resent time 
by introducing additional simj)lifications. One method 
proceeds by aj)proximating uw with some simjile functions 
having C(‘rtain jiarameters temporarily unspecified. Values 
for tlie hitter are determined by substituting the functions 
into the integral ecpiation and satisfying the ecpiation at a 
limited number of ])oints ecpial to the number of unspecified 
parameters. At this point in the analysis, Oswatitsch 
assumed that the vai-iation of uw across the chord could 
be represented by a parabola, one or two lialf-parabolas, 
or a rectangle combined with a half-parabola, as illustrated 
in figure 15, all of unspecified height and chord wise extent. 



ITci'RE 15. — Assumod chordwise variations of iiw used by Oswatitsch. 


Xo account was taken in the integral of the iidluence of 
the region between the leading edge and the station where 
o:i (he forward ])art of the airfoil nor between the 
station wluu’e UL^y=^) on the rear of the airfoil and the 
trailing edge. Uj)on inserting a selected combination of 
the abovi' mentioned elements into eciuation (6()) and inte- 
grating, there resulted a system of simultaneous quadratic 
algebraic equations having as many members as there 
were elements in the selected general form of solution. 
In many cases, Oswatitsch used only one element, either 
a ])arabo]a or a half-])arabola and assumed a mean value 
for h for tlie entire chord. The method included no provision 
for the imj)j’ovement of the result through iteration or otlier 
ni(‘ans, the only measure of the accuracy being the degree 
of corresj'ondence between the initial and the final Uw dis- 
tribution^=. Nevertheless, the results presented in references 
20 and 21 were encouraging in that they showed many of 
the known characteristics of transonic flow about airfoils. 
In ])articular, shock waves apjieared when the critical 
Mach numl)er was exceeded and moved rearward with 
further increases in Mach number. However, because tl\e 
initial valiu's for Uw wiu*e generally sul)stituted not only 
into the integral but also into the term Uw^j2 standing 
outside tl e integral, and because the tangency iTCpiin'inent 
on the I and L curves necessary for smooth transition 
from subsonic to suj)ersonic velocities was not realized, a 
multiplicity of solutions was often obtained. In one case, 
three soli tions were actually given and more could have 


been obtained which would have fulfilled equally well the 
condition of corresjiondence between initial and final result. 
In addition, the quality of the results appeared to deteriorate 
when the Mach number was increased to higher values, the 
uj)j)er limit of acceptability aj)j)earing to be a Mach number 
of about 0.88 for 6-percent-thick airfoils. 

The integral equation method has been developed further 
in references 22 through 25 by Gullstrand. In the first of 
these, reference 22, e([uation (61) is simplified to a single 
integral ecjuation through the use of a more elaborate veloc- 
ity-distribution function tlian that given in equation (63), 
and the resulting e([uation is solved by an iteration process. 
The introduction of an iteration procedure makes a marked 
improvement over the metliod of Oswatitsch since it then 
becomes ])ractical to increase greatly the number of elements 
with which is represented. The method is applied to 
three 6 percent-thick XACA low-drag airfoils and the re- 
sulting velocity (listril)utions are given. In common with 
the original method of Oswatitsch, difficulties occur when the 
Mach number is too close to unity. The highest Mach 
number for which results are given is 0.91. More recently, 
Gullstrand has ])resented approximate solutions for the 
velocity distriliutions on symmetrical airfoils in sonic flow 
in reference 23.^ In both of these works, however, the itera- 
tion process proceeds by inserting the known values into 
both the integral and the term Uw^j2 standing outside the 
integral. This procedure is equivalent to replacing the 
second-degree equation for u^r with a linear equation and 
obscures or loses many of the characteristics of the quadratic 
solution discussed in the jireceding section. 

The ])resent analysis also proceeds through the use of an 
iteration scheme, ])artly numerical and partly graphical, but 
the known values are substituted only into the integral at 
each stej) of the process. The quadratic nature of the equa- 
tion is thereb}^ retained and the tangency condition on the 
1 and L curves can lie fulfilled. Satisfaction of the latter 
requirement is essential for uniqueness and convergence. 
Inasmuch as the calculations for the circular-arc section 
were well advanced when Gullstrand’s pajiers were received 
and were producing reasonable results witli the simple 
velocity-distribution function of Oswatitsch, it was decided 
to continue rather than to start over using a more elaborate 
velocity-distribution function. It has lieen found, however, 
that additional attention must be paid to the influence of 
the region behind the trailing edge at the higher Mach num- 
bers. l"pon observance of these additional refinements, re- 
sults are obtained for all Mach numbers up to unity. At the 
lower Mach numiiers, these results are in general agreement 
with those found by the simpler, alt hough more approximate, 
method of Oswatitscli. The present method carries right 
on, however, into the higher Mach number range where the 
simpler method met with difficulties and succeeds in showing 
the well-known invariance of local Mach number with 
changes in the free-stream Mach number as the latter 
approaches unity. 

2 Since eqmition (61) degenoratos at a Mach nuinber of unity where /3=0, the sonic results 
of reference 23 are obtained by first developing a new double-integral equation to replace 
equation (61) and then simplifying and solving by an iteration procedure. The integral 
equation is only used, however, to determine the velocity distribution over the portion of the 
airfoil forward of the station of maximum thickness. The remaiuder of the solution is ob* 
tained by means of the theory of characteristics. 
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In tho present calculations for circular-arc airfoils, the 
necessity for additional refinement l)egins at Mach niimhers 
somewhat greater than that at which the shock wave reaches 
the trailing edge. Consequently, no attempt is made to 
account for the influence of the region behind the wing until 
the Macli number becomes sufficiently large for the sliock 
wave to reach tlie trailing edge. At higher Mach numbers, 
the influence of this region is api)roximated in the following 
manner. First, it is assumed that the shock wave which 
stands at the trailing edge remains of the strong family as 
tlie free-stream Mach number is increased to unity. It is 
assumed furthermore that the flow is parallel to the x axis at 
the shock position and that the shock wave is normal to the 
local flow. With these assumptions, it follows from equa- 
tion (45) that the values of u immediately behind the shock 
wave are related to tliose immediately ahead of the shock by 

Ut = 2~Ua ( 68 ) 


Since Ua is given by equation (63), it follows that Ui, can be 
expressed in terms of the values of Uw «ud b iinmediatel}" 
ahead of the shock wave in the following manner: 


u,{z)-2 + 


(09) 


Because the shock terminates within a finite distance of the 
airfoil at all subsonic Mach numbers, equation (69) is 
appropriate only within a certain range of |?|, namely, 
that for which ^a>l- In the present calculations, the end of 
the shock occurs at such large [z\ when the shock is at the 
trailing edge that equation (69) has been used to represent 
the values for u behind the trailing edge for all [z\. It is 
further assumed that the contribution of the region behind 
the trailing edge can be satisfactorily approximated by 
equating u to ud for all points behind the trailing edge. It is 
recognized that neither of these assumptions constitutes a 
good approximation for u at great distances from the wing, 
but the attenuating influence of distance diminishes the 
error in the values of the integral at points on the airfoil 
surface. In this way, ec{uation (66) comes to be replaced 
with the following relation determined from equation (61) 
by performing an integration with respect to J under the 
assumption that the variation of u witli ~z (or 1) is given by 
equation (63) for stations ahead of the trailing edge and by 
equation (69) for stations beliind the trailing edge. 




where E is as defined by equation (67) and F is given by 



4 

TTd+X^^ 


[x|X|(3-X^- 


2(l-3.Y2)lnlA^l~(3+2X- 


-A^^)] (71) 


The nature of F{X) is illustrated graphically in figure 16. 
In this case it is apparent that the three integrals of equation 
(70) taken togetlier corresiiond to 7/2 in equation (56). 

To summarize, equation (06) is used in the present calcu- 
lations when the Mach number is less than that at which the 
shock wave first reaches the trailing edge, and equation (70) 
is used for higher Mach numbers. No significant discon- 
tinuities are jirodiiced in the value of the integrals, however, 
since the contribution of the additional integrals only 
becomes significant at Mach numbers greater than (hat at 
which the jirocedure is changed. 



NUxMERICAL EVALUATION OF INTEGRAL 


One of the ])rincipal ste])s in the iteration method used 
herein for the solution of equations (66) and (70) is the 
evaluation of the integrals. Since uw and h are generally 
prescribed by a set of numerical values rather than by 
analytical functions, a numerical technique has been used 
for the integration. This process consists of replacing the 
prescribed uw distrilmtion with a stepwise approximation as 
indicated in figure 17, introducing a mean value for h for 
each of the rectangular elements, integrating to determine 
the contribution of a single element, and summing the 
influence of all the elements. The contribution of a single 
element of width I situated on the airfoil chord, as typified by 
the shaded area of figure 17, is given by 



( 72 ) 
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Figure 17.— Stepwise approximation used in evaluation of integral 
''vhen shock wave is forward of trailing edge. 


l"pon perfo •manco of the indicatod o])orations, the following 
expression is obtained for/i: 

1 2A In 1.4| -.4(1 +.4') [(1 +AY+12] + 

1 fSw B 


2 \B\ 


[(l+/l^/-(l+/l')'+8(l + /yV8] + 


12/^(/;2-l)ln|/i|-«(l+/i2)[(l+/i2)2+i 


4 


where 


1:A2{x-Q \ (2h\ I (-x-iA/2k\ 

2hV “4V6J^2V h )\hi) 

I, -2 (i-li) {^\ 


B= 


2b, 


li / \ / 


(73) 


(74) 


Thus, the integrals in equations (66) and (70) that eontain 
E aj’c approximated as follows: 


/ % f 6') [( V')' (tl 


(75) 


It is evichiiit from equation (70) that the eontribution of 
a single element of width I situated behind a shoek at the 
trailing edg(', as typified by the shaded area of figure 18, 
is composed of two parts. The first depends on an integral 
involving E{X) and is evaluated using /i in the manner just 
described for elements on the airfoil. The second deiiends 
on an integral involving F{X) 


gi 




d\ (7(1) 



Figure 18. — Stepwise approximation used in evaluation of integral 
when shock wave is at trailing edge. 


which, upon evaluation, yields 

4 [(l+^l^^+(l+-4^ — 2]— 2Alln \A\ — 

(77) 

where A and B remain as defined in equation (74). Thus, 
the integral which contains F in equation (70) is approxi- 
mated as follows: 


27^11 1 /?|-/^ ( 1 + 7 ^ 2 ^ 


} 




(78) 


Values of /i and (ji are presented graphically in figures 19 
and 20. (It is noted that the graphs in reference 21 that 
correspond to figure 19 of this report arc mislabeled.) 

With the simplifications introduced in this section, the 
function / of equation (56) is approximated by 


I(x)=2^ 




(79) 


for Mach numbers less than those for which the shock wave 
is situated at the trailing edge of the aii’foil, and 


I{J) = 


2{2=jv,[(vi-)(i); 

“i--sy{(^)(i),(|)]_ 




(80) 


for larg(?r Macli numlicrs. 
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Figure 19. — Variation of influence function /i with 

Computed from equation (78). 


-r- for various values of 


DETERMINATION OF 

The term Ul^ that appears in the integTal equation for 
transonic flow represents the values for u given hy linear 
theory for points on the airfoil surface. Its values can be 
obtained from the general two-dimensional solution for Ul 
given in equation (47) by performing the indicated operations 
and setting z=(). As noted in equation (51), simplification 
occurs for the symmeti*ical nonlifting airfoils being considered 
herein because 

Au=Uu^Ui=0 


^ d c>r ^ (le ^ 


e>f 


df 


(81) 


tlius 

Ul^ 






(82) 


The Cauchy principal value is understood in the last integral. 
The present calculations are for thin circular-arc airfoils for 



Figure 20. — Variation of influence function < 7 i [('T"'KS] with 
2/ • X — ^ i 

for various values of y— • Computed from equation (77). 

which the i-ednced ordinates Z„ of the u])i)ei' surface are given 
hv 




(83) 


where r represents a reduced thickness i*atio which is related 
to the actual thickness ratio as follows: 


-_kUo t 
^ 3^ c 


(84) 


Performing the indicated operations gives the following 
expression for UL^r' 




(85) 


We thus have, l)v sidistitnting equation (85) into equation 
(59), 


UJ'=25v-1=!'[i+(H)1"o-j]-‘ 


( 86 ) 


It is seen that r plays the role of a similarity parameter. 
Thus, for a family of symmetrical nonlifting airfoils hav- 
ing the same thickness distribution depends only on r and 
the position coordinate 7/c. Inspection of the integral equa- 
tions for transonic flow shows that their solutions for also 
depend only on 7 and 7/c. 

Many previous pajiers on transonic flow (e. g., refs. 12, 
14, 16, 17, and 18) have used tlie symbol to designate a 
Mach number thickness ratio parameter diflerent from the 
parameter 7 used herein. The definition of suitably gen- 
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cralizod to allow for various exprossions for k, and tlio rcda- 
tion l)(‘tw(‘( n <^*id r are as follows: 


^o=- 


] 



(87) 


The ])arainot('r ])ossossos tllo distinct practical advantag-c 
of ap|)roaching zero rather than infinity as the Mach nuniher 
aj)])roach('s unity, ddie corresponding ])arameter refei*ring 
to the local conditions has coniinonly been designated by 
the symbol 




\-\p 

[UMrnY'^ 


( 88 ) 


where M is the local Alach number. Since the (luantity 
1— d/“ is siiown in ecpiation (24) to be ecpiivalent, in tlie 
])resent aj)pi*oximation, to \—Mo^—ku, the redation for ^ 
may be rewritten as 


\-Mo^-ku 

[Uokmf^ 


(89) 


which may be expressed in terms of u and r as follows: 




u—\ 

72/3 


(90) 


ITERATION SOLUTION OF INTEGRAL EQUATION 

Solutions of the simplified integral equations have been 
obtained for transonic flow about circular-arc airfoils by 
using a numerical and grajihical ])rocess. Four slightly 
different technicpies are used depending on tlie Mach num- 
ber or, mo ‘e precisely, on the value of r. One jirocedure 
is used for th(‘ subcritical Mach number range, a second is 
used in an atlemjit to find sujiercritical sliock-free solutions, 
a tliird to determine supercritical solutions in the range of 
T for whicli the shock jiosition is forward of the trailing edge, 
and a fourili for still larger values of r. The latter range 
(‘xtends uj) to a free-stream Mach number of unity. 

Subcritical flows. — Solutions for small Mach numbers can 
be obtained by a direct iteration process starting with the 
linear-theory solution. In detail, the calculations proceed in 
the following manner. The values of corresponding to a 
given T (s]»ecified thickness ratio and Mach number) are 
calculated from equation (85). The curve is apjiroxi- 
mated witl a stepwise distribution of ten stejis, and the 
values of the I curve are computed therefrom using equation 
(79). Knowing the values of tTz.^^.(x) and 7(x),one obtains a 
first approximation to i^ir(^) using ecpiation (59). The 
jirocess is now rejieated using the values for uir to calculate a 
new I curv(‘, from which a second ajiproximation can be de- 
termined forT'H., Ptv. A typical set of results illustrating the 
convergence of the jirocess is shown in figure 21. The 
process apjK'ars to converge rajiidly, at least in the piTsent 
calculations, whenever the Mach number (or more properly 7) 
is sufficient Iv small that u does not exceed unity (or the local 
v(‘locity does not IxTonie supei'sonic) at any jioint in the 
calculation. If 7 is increased to a value such thatTZ becomes 
ecjual to one at any stej) in tlu' iteration, however, the jirocess 
terminates abru])tly with the ap])earance of complex values 
of Ti. Sinc(‘ the results of successive iterations oscillate in 


this Mach number range and since the starting distribution 
])rovides maximum values that are too small, any such 
termination occurs in the first iteration stej). Such an 
abrupt termination of the calculations is in marked contrast 
to the well-known jiroperty of the classical iteration methods 
of ])roducing second- and higher-order solutions indicating 
shock-free mixed flows. It will be shown in the following 
discussion that the difference in behavior is not a product of 
the additional assumjitions and apjiroximations introduced 
herein, but stems directly from the recognition and retention 
of the (piadratic nature of equation (bl). 



' 0 .2 .4 .6 .8 1.0 

x/c 


Fkjure 21. — Results of iteration calculations in the subcritical range, 

7 = 0.229. 


This point can be discussed in greater detail by confining 
attention to only the first step of the iteration procedure in 
which the / curve is calculated using Uw=UL^y If the present 
iteration procedure is used, if the UL^r curve is replaced with 
ten rectangular steps as shown in figure 17, and if the / curve 
is calculated using equation (79), then equation (66) yields 
the following result for at the 50-percent-cliord station 
of a circular-arc airfoil: 

(I ' i -y/i 1 (91) 

For low Mach numbers, 7 is small, the discriminant is posi- 
tive, and the first approximation for u^r can be readil}' ob- 
tained using the minus sign. As 7 increases, the value of 
the dis(*riniinant decreases and becomes negative wdien 7 
exceeds 0.490. This value of 7, therefore, rejiresents an 
ujiper limit for the Mach number at which the present 
])rocess can lead to a useful result if the starting Uiy distribu- 
tion is taken to be that given by linear theory. 4'he value 
of Uw{cl2) given by eciuation (91) is equal to (47 /tt) at 
small 7, and attains a maximum value of unity when 7=0.490. 

The classical iteration methods of subsonic flow theory 
have not shown such an abrupt termination of solutions 
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upon attaininent of sonic velocity. The same situation 
occurs with the present ecpiations if the iteration procedure 
is altered so as to conform with the classical methods. This 
change consists simph' of substituting ULy^, into not only tlie 
integral of equation (66) but also the term standing 

outside the integral, thereby converting the quadratic 
equation for Uw into a linear equation. If this procedure is 
adopted together with the same value for the integral that 
was used to obtain eciuation (91), the result of the first 
iteration is 

«.^(0=^^+«-295T^ (92) 

Since a linear equation is solved at every step of the iteration 
process, the procedure never terminates. Values for Uw in 
the midchord region become larger with every iteration step, 
however, and it seems to be the prevailing belief that con- 
vergence is obtained only for Mach numbers less than the 
critical. 

It is of interest to compare the approximate results of the 
first iteration step reviewed above with the exact results for 
the same airfoil given by Hantzsche and Wendt in reference 
32. If the latter results are made comparable to the present 
results by taking the limiting value corresponding to small- 
disturbance transonic flow theory, the exact result of the first 
iteration step is 

( 1 )=^ ^ 

The first term of both equations (92) and (93) is that given 
by linear theory and is the same in both calculations. The 
difference in the coefficient of the second term is the result 
of the errors introduced in the approximate solution of 
equation (61) (i. e., the velocity-distribution function, finite 
steps for Ti-wj etc.) and is some sort of a measure of the accu- 
racy of the approximate calculations. Just as with equation 
(92), a value for Uw can be calculated for all t, although the 
question remains of whether the result is a valid first step in 
a convergent process. 

If Hantzsche and Wendt liad performed their iteration 
calculations in a manner comparable to that described herein 
so that the values of Ul iht not introduced as an approxima- 
tion for u in determining the important influence of the points 
lying near P, they would have obtained the following relation 
(again to the approximation of transonic flow theory) as the 
result of the first iteration: 



This result is directly comparable with that given in equa- 
tion (91) except that the double integral of equation (61) is 
evaluated exactly rather than approximately. Again, the 
accuracy of the approximations can be evaluated by com- 
paring the two equations. More important, however, is the 
fact that the exact results also terminate when r exceeds a 
certain value (0.503) and that the corresponding value for u 


is unity. This subject is discussed in greater generality in 
Appendix B. 

As noted above, the first step of the ])resent iteration pro- 
cedujT starting with Uj^ jirodiices maximum values for u 
that are too large. Consequently, the calculations terminate 
at a value of r that is somewhat less than corresponds to 
the true critical Mach number. vSolutions can be obtained 
for the remainder of the subcritical Mach number range, 
however, if the starting Uw is selected having the same general 
form but larger values than the corresponding Ui^^y distil- 
bution. A good starting distribution can be obtained by 
simply extrapolating the final results for smaller r. 

Calculations of tlie type just desci’ibed have been cariled 
through for several values of r less than the critical value 
(7^^= 0.598), defined as the smallest value of 7 at which sonic 
velocity {y^r=l) is attained in the final transonic solution. 
The results of these calculations are presented in the form of 
chordwise 17^' distributions for various 7 in figure 22, and in 



Ficjure 22. — Chordwise variation of for various r in the subcritical 

range. 

the form of curves of versus 7 for various xlc in figure 23.^ 
The corresponding values given by linear theory are also 
shown in figure 23. Comparison of the two sets of curves 
reveals that the values of Uw eoincide witli those for 

3 The results for t = tct are also included on these graphs in order to complete the subcritictil 
range, even though they are obtained using a diilerent iteration procedure. Since the pro- 
cedure is the same as that described in the next section, however, no further comment is 
necessary at this point. 
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Figure 23. — Variation of uw and with f for various values of xic. Subcritical range. 
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small r (low Afach number) but become larger than those for 
at most stations on the airfoil as r approaches Tct- "riiese 
trends, which are in accord with generally accepted experi- 
mental and theoretical results, will be discussed at greater 
length in a later section. 

Shock-free supercritical flows. — Although it is indicated 
in the preceding discussion that the present iteration process 
does not lead to supercritical shock-free flow about airfoils 
when the starting point for the iteration calculation is Ul^ it 
is considered of importance and interest to ascertain the be- 
havior of the present approximation when such a solution is 
deliberately sought b}^ starting with sufficient^ large values 
for u that I is greater than L. The interest in this matter is 
heightened by the fact that the less detailed method used by 
Oswatitsch leads to such results over a limited range of super- 
critical Mach numbers. (See figs. 4 (c) and 4 (d) of ref. 20 or 
figs. 7 (c) and 7 (d) of ref. 21.) Accordingly, iteration calcula- 
tions have been performed starting with r greater than 
and a symmetrical shock-free u^r distribution containing a 
region of supersonic flow (^^^>1) over the middle portion of 
the chord. A slight change in the iteration procedure is 
necessary, however, to eliminate the difficulty arising from 
the requirement described in. an earlier section that I=L 
at the sonic points. To carry out the solution in the same 
manner as before would require that the distribution be 
found for which the corresponding I curve is tangent to a 
given L curve. Kather than attempting to find solutions by 


such indirect means, a method is used in which a new value 
of T is determined in such a manner that the tangency con 
dition is fulfilled. In particular, the procedure is to calculate 
the / curve using equation (79) and the assumed values of r 
and Uiy. The next step consists of plotting the I curve and 
fitting an L curve computed from equation (86) for whatever 
value of T is necessary to fulfill the tangency condition as 
shown in figure 10. A new set of values for Uy. can now be 
calculated using equation (59). 

U\y=lzL->Jl (69) 

where the minus sign is used at stations upstream from the 
forward sonic point and downstream from the rear sonic 
point and the plus sign is used for the intermediate stations. 

This process has been carried out for a number of assumed 
initial Uw distributions. Figure 24 shows a typical set of 
results obtained by starting with an initial value for r of 0.65 
and a Uw distribution obtained by extrapolating the trends 
indicated by the solutions for subcritical flows. The as- 
sumed initial distribution is indicated by a dotted line in the 
upper portion of part (a) and the associated / and L curves 
are shown immediately below. The tangency requirement 
is fulfilled by taking r = 0.642 which compares well with the 
initial value of 0.65. The corresponding Uw distribution 
calculated therefrom is shown in the upper portion of part (a). 
Several points of interest are to be noted. The first, of 
course, is that tlie assumed and calculated Uw distributions 



x/c 

Figure 24. — ^Uesults of successive iteration steps in an attempt to find a shock-free supercritical solution. 
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are distinctly dissimilar. The seco/id is that the calculated 
distribution bears a marked reseml)lance to a residt given 
by Osvvatitsch (see hg. 4(d) of ref. 20 or fig. 7(d) of ref. 21) 
which may be described as a s^unmetrical Uw distribution 
containing an abrupt ev'pansion on the front of tlie airfoil 
and a compression shock symmetrically situated on the rear 
of the airfoil. 

Since the assumed and calculated Uw distributions are too 
dissimilar to be regarded as solutions, the itei’ation process 
was continued, using the calculated results of the first step as 
the initial values for the second step. The results of this 
calculation are summarized in part (b). As indicated, r 
changes from 0.642 to 0.671 and the new Uw distribution is 
again markedly dissimilar to the initial distribution. The 
calculations have been carried through several more steps 
and the results are shown in parts (c) through (f). It can 
be seen that the initial and calculated Uw distributions do not 
agree, even qualitatively, until several steps of the iteration 
process have been completed. In the meanwhile, the value 
of r has decreased to approximately that for the critical 
Mach number and the region of supersonic flow has effec- 
tively disappeared. The sixth and last step of the iteration 
calculations presented in figure 24 has produced a value of 
r = 0.600 and a Uw distribution nearly identical with that 
shown in figure 22 for the critical Mach number (rcr= 0.598). 
Although this figure shows the results of only one series of 
calculations, similar results have been obtained starting with 
other assumed symmetrical, shock-free, supercritical Uw dis- 
tributions. No case was found in which the calculated values 
repeated the assumed values until r had decreased to ap- 
proximately Tct and the sujjersonic zone had vanished. 

It should be noted before leaving this section that these 
results are not presented with the intention of proving or 
disproving anything about the more fundamental questiozi 
of the existence or nonexistence of shock-free transonic flows. 
The purpose, i-ather, is to illusti'ate the behavior of the 
present approximation furnished by the simplified integral 
equation and the iteration method of solution. 

Supercritical flows — shock wave forward of the trailing 
edge. — ^The preceding section has shown how the present 
method of calculating velocity distributions on thin air- 
foils fails to converge when 7 is greater than ^ind the flow 
is assumed to be shock free. The identical iteration process 
will converge rapidly to a solution, however, if the initial uw 
distribution contains a discontinuity in accord with the 
shock relations. The method starts by selecting a value for 
r and assuming a reasonable distribution for Uw over the 
chord. The main point to observe in the selection of uw is 
to include a shock wave (a discontinuity in uw) through which 
jumps from 1+A immediately ahead of the shock to 1— A 
immediate^ behind the shock. As noted previously, such 
a jump in uw is consistent with the assumption that the 
shock wave is a weak noi*mal shock. The next step is to 
calculate the I curve using equation (79) and the assumed 
values of and 7. The / curve is plotted and an L curve 
is computed using equation (86) and whatever value of 7 is 
necessary to fulfill the tangency condition, as illustrated in 
figure 13. A new set of values of uw can now be calculated 
using equation (59), taking propei* care to change from the 


minus to the plus sign at the point of tangency and then re- 
turn to the minus sign aft of the assumed shock position. In 
this way, a new approximation for tlie Uw distribution cor- 
responding to the new value for 7 is obtained, but the posi- 
tion of the shock wave is unaltered. If the new values for 7 
and Uw are sufliciently close to the assumed values, it is pre- 
sumed that an approximate solution has been found. In 
general, however, such a close correspondence is not obtained 
after the first step, and the entire calculation is repeated 
using the new values for 7 and uw in place of those assumed 
initially. Throughout the process, the position of the shock 
wave is kept fixed, and the value of 7 is allowed to vary as 
necessary. In this way, the iteration process selects the 
solution for a given shock position rather than for a given 
Mach number. This process was carried out with the shock 
wave situated at several diflerent stations on the airfoil 
chord and was found to converge rapidly even when the 
selected initial values for uw aiid 7 dift'ered considerably from 
their final values. A typical set of results (namely, that for 
the shock position fixed at 90-percent-chord station) is shown 
in figure 25. In common with the rest of the calculations, 
tiie initial values for uw and 7 were selected by extrapolating 


Initial (t= 1.043) 

fFirst approx, (7=1.033) 
[Second approx. (f= 1.029) 
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Figure 25. — Resiilt.s of iteration calculations with shock wave at 90- 

percent chord. 
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the trends indicated hy the calculations for more forward 
sJiock positions. 

Calculations of the type just described have been carried 
through with the position of the shock wave fixed at 60-, 
70-, 80-, 90-, and 100-percent chord. The results are 
presented as chordwise u^r distributions for various r in 
figure 26 and as the variation of Uyy with 7 at selected values 
of 3r/c in figure 27. 



x/c 

Figure 26 . — Chordwise variation of for various f, shock wav'e for- 
ward of trailing edge. 


In order to test the method further, calculations were 
repeated in many of the cases using a difi’erent set of initial 
values for 7 and u^r. The results of two such calculations 
for the shock position at the 90-percent-chord station are 
shown in figure 28. Although the results sliown in part (a) 
were obtained using values for uw that were purposely 
selected to be too large and those of part (b) to be too small, 
it can be seen that similar results are obtained after only a 
few steps of the iteration process. These same results are 
presented in a second manner in figure 29 wherein the values 
of u^r at several stations on the chord are plotted as a func- 
tion of 7. The final values to which the calculations con- 
verge (i. e., those given in fig. 27) are indicated by the solid 
line. The points connected by the dotted lines are the 
values obtained at eacli step of the calculations. It may 
be seen that the present procedure appears first to place 
the values of and 7 on the curve of correct solutions and, 
subsequently, to converge to the final result. 

Supercritical flows — shock wave at the trailing edge. — The 
calculations just described indicate that the shock wave 
moves rearward with increasing r and reaches the trailing 
edge when 7=1.34, corresponding to a Mach number of 0.92 
for a 4-percent-thick circular-arc section. The calculations 
for larger subsonic Mach numbers were performed under 
the assumption that a strong shock wave remains at the 
trailing edge. As in the preceding analysis, the shock wave 
is assumed to stand perpendicular to the x axis, and the flow 
is assumed to be normal to the shock wave, even though 
these conditions cannot be correct at the base of the shock 
wave since the flow must turn through a finite angle. 

The method of calculation used for the higher Mach 
number range where the shock wave stands at the trailing 
edge is essentiall}' the same as that used for mixed flows at 
lower Mach numbers. A slight modification enters in that 
the / curve is calculated by means of equation (80) rather 
than (79). This change is made because the contribution 
of the region immediately behind the trailing edge lieconies 
of increasing significance as the Mach number approaciies 
unity. Since tlie contribution is small at the lowest Mach 
number for which the shock wave is situated at the trailing 
edge, no significant discontinuity is produced in the results 
by this change in procedure. A difficulty arises in the itera- 
tion process, however, because there is no longer any dis- 
tinguishing feature to fix the value of 7 for which tlie solution 
is being sought. To review, for pure subsonic flows, 7 itself 
can be maintained at a fixed value from step to step in the 
iteration process. At Mach numbers somewhat greater 
than tlie critical, where the shock wave stands on the airfoil, 
the tangency requirement makes it difficult to prescribe 7 
directly but a satisfactoiy method is obtained by fixing the 
sliock position and carrying out the iteration process until 
the associated value for 7 is found. At still larger Mach 
numbers, the first method is of no avail and the second 
method cannot produce a unique result since the shock wave 
is considered to be fixed at the trailing edge over the rest of 
the range of subsonic free-stream Mach numbers. For lack 
of a better method, the calculations were carried out for 
this range using a procedure that might be described as 
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Approx r 



(a) 


0 .2 .4 _ .6 .8 1.0 

x/c 

(a) Initial uw too large. 


Approx T 

Initial l.lOO 

First .864 

Second .975 

Third .974 

Fourth .997 

Fifth 1.006 


Sixth 1.0 15 



(b) 


0 .2 .4 .6 .8 1.0 

x/c 

(b) Initial uw too small. 


Figure 28. — ^Resnlts of iteration calculations started with deliberately poor assinnption for , shock wave at 90-percent chord. 


being partly iteration and partly trial and error. The 
calculations proceed as follows: A value of 7 is selected 
and an estimate made for the associated Uw distribution. 
On the basis of these values, calculations are performed 
resulting in a new set of values for 7 and Uw Similar cal- 
culations are repeated on a trial-and-error basis using 
the same Uw distribution but various values for 7 until 
the resulting 7 is equal to the assumed value. These 
values are then plotted on a graph of Uw versus 7. If these 
values were an actual solution of the integral equation, the 
resulting values of Uw would also be equal to the assumed 
values. In the present calculations, however, the resulting 
values of uw are generally found to be somewhat smaller 


than the assumed values. These new values for Uw, to- 
gether with a smaller value for 7, are next taken as the start- 
ing values for a second series of trial-and-error calculations. 
Again, the Uw distribution is held fixed as various values are 
tried for 7. The calculations are again terminated when 
the value of 7 is so selected that it repeats itself. The values 
of Uw and 7 are then plotted on the graph. Unfortunately, 
it was not found possible in the present calculations to de- 
termine a set of values for uw and 7 that would repeat them- 
selves precisely. In all cases investigated, the values of 
Uw and 7 were found to diminish somewhat in successive 
iteration steps. Tiie source of this difficulty has not been 
ascertained at the present time. It could be due to one or 
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.6 .8 1.0 1.2 


Fkjure 29. — Results of calculations with shock wave at 9()-percent 
cl ord. (Cross plot of data in figures 25 and 28.) 

more of tl e simplifying assumptions introdiieod in the ap- 
proximate calculations or it could be due to the fact that 
an iteration process was not devised which would lead to a 
definite result, as at smaller 7. 

There a:’e a number of points, however, which tend to 
indicate that the successive values of uw hiid 7 obtained 
after the lirst few iteration ste])s may l)e regarded as near 
solutions. One of these concerns the fact that calculations 
made starting with different initial values for uw and r 
converge to a common result after the first few iteration 
steps. To illustrate, results of typical series of calculations 
starting with 7 of about 20 and three different uw distribu- 
tions are sliown in figure 30. It can be seen that the three 
sets of calculations all converge to determine a single line 
after the first few iterations. In order to promote insight 
into the significance of the indicated variation of 7 from ste]) 
to ste]) of the iteration process, the corresponding Mach 
numbers or a 4-percent-lhick section are also shown. 
Comparison on this basis shows that the variation of Mach 
number frDin step to stej) of the iteration process is very 
small. 



Fkjure 30.— Results of iteration calculations started with deliberately 
poor assumption for u\v, shock wav^e at trailing edge. 

An additional factor lending credence to the ap|)licability 
of the solutions for large 7 is provided by the phenomenon 
of the Mach number freeze wherein the local Mach number 
is invariant with changes in the free-stream Mach number 
when the latter is near unity or, more j)recisely. 



It is known from the ])a])ers of Vincent i and Wagoner (ref. 
14), Liepmann and Bryson (refs. 28 and 29), and others that 
the corresj)onding ajiproximate relation yielded by the small- 
(listurl)ance transonic theory is 



If the parameter k involved in the definitions of J is inde- 
pendent of as is the case when k is taken to be as given 
in equation (17), and if the local Mach number M is cal- 
culated using eciuation (24), the above two relations are 
com])letely e(|ui valent. If k contains d/^, however, as it 
does in the preferred definition given in equation (7), (L\fldM^ 
only vanishes when d/, as well as d/^, equals unity. If it is 
assumed that the freeze extends over a finite range of Mach 
numbers, the variation of u with 7 must obey the following 
relation for large 7: 

^=:^^=const. (97) 
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Equation (97) has been used to compute a curve of u^y 
versus r under tJie assumption that the freeze extends at 
least to as small a value of 7 as 10, which corresponds to a 
Mach number of 0.978 for a 4-percent-thick section. These 
lines, which all have the characteristic form 


Uw—l 

{uw)t=\q 1 \ ib / 


(98) 


are also shown in figure 30. It can be seen that the lines so 
determined are almost identical with those determined previ- 
ously after the first few steps of the iteration process. This 
comparison shows that the solutions obtained by the iteration 
process possess the phenomenon of the Mach number freeze. 
A further check was made by extrapolating the curves of u^y 
versus 7 to 7=40 with the aid of equation (98) and using 
the values thus obtained as the initial values for an iteration 
calculation. It is found that these values nearly satisfy the 
integral equation although the values of u^r and 7 diminish 
slightly from step to step during the iteration ])rocess, just 
as at the smaller values of 7. The significance, in terms of 
Mach number, of these changes is even smaller than at 
lower Mach numbers. For example, the changes produced 
in the values of u and 7 in the first step of the iteration 
process correspond, for a 4-percent-thick section, to an 
indiscernible change in the local Mach number and a change 
in the free-stream Mach number of less than 0.001. On the 
basis of these results, it appears proper to extrapolate the 
curves of versus 7 to values of 7 approaching infinity, 
corresponding to a Mach number of unity, by using equa- 
tion (98). 

The results of the calculations for the range of 7 for wliich 
the shock wave is at the trailing edge are summarized in 
figures 31 and 32. 

RESULTS 

RESULTS IN TERMS OF REDUCED QUANTITIES 

The calculations described in the preceding sections have 
produced values for the velocity distributions at the surface 
of tliin circular-arc airfoils in flows having free-stream Mach 
numbers ranging from zero to unity. These results are pre- 
sented in graphical form in figures 22, 23, 26, 27, 31 and 32 
in terms of the reduced quantities u and 7 defined in equations 
(29) and (84) and repeated below. 

k - Uok t 


where k represents the coefficient of tlie nonlinear term in 
the approximate differential equation for (p (eq. (6)). 

The results of a number of previous investigations of 
transonic-flow theory have been presented in terms of the 
quantities ? defined in equations (87) and (88) 

, l-M^ 

[Uokmr^' ^ 

The relations between the two sets of quantities are given 
by equations (87) and (90) 




x/c 


Figure 31. — Chordwise variation of uw for various t, shock wave at 

trailing edge. 


Pressure distribution. — In many applications, the quantit}^ 
of prime interest is the pressure distribution rather than the 
velocity or Matdi number distribution. Because of the 
simple relation between the pressure coefficient and the 
perturbation velocity provided by equation (16), however, 
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(a) f<4.o 

Figure 32.— Variation of Uw with ? for various values of x/c, shock wave at trailing edge. 
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(b) Large r. 

Figure 32. — Concluded. 
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it is a simple', matter to determine the pressure distribution 
once the velocity distribution is known; thus, 


2u_ - 2{tlcyi^ ( u \ 

Uo Uok qJoky'Ar^'V 


(09) 


The latter expression for Cp suggests the introduction of a 
reduced pressure coefficient (\ defined by 


-_(Uoky'^ 



( 100 ) 


Equation (100) may he rewritten in terms of ? as 

follows: 

C,= -2{i-U ( 101 ) 

The foregoing results have been used to calculate Cp for 
numerous stations on the airfoil chord for various (or 7). 


The results so obtained are presented in figure 33 in the form 
of chordwise pressure distributions {Cp vs. xjc) for various 
^ 0 . This figure is presented in three parts: Part (a) contains 
the results for subcritical Mach numbers, part (b) for slightly 
supercritical Mach numbers for which the shock wave stands 
on the airfoil surface, and part (c) for still larger Mach num- 
])ers for which the shock wave is at the trailing edge. 

The same results are presented in another manner in figure 
34 in which the variation of (\ with is plotted for various 
stations xjc along the airfoil chord. This form of presenta- 
tion is the counterpart, in terms of reduced quantities, of the 
plots commonly found in many experimental investigations 
showing the variation of pressure coefficient Cp with free- 
strcain Mach numl)er Mq at prescribed points on the surface 
of wings and bodies. The corresponding curves given by 




Figure 33. — Chordwise variation of Cv ^r various ^o. 


Figure 33. — Continued. 
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linear theory are shown in figure 34 by dotted lines. Tliese 
curves are computed using the equation 

obtained by direct substitution of equations (87) and (100) 
into equation (85). It can be seen that the present results 
and those of linear theory are in good agreement for values 
of ^0 considerably less than zero. For Jo near to zero, how- 
ever, the present results display a l)ehavior completely dif- 
ferent from that indicated by linear theory. This is as it 
should be since it is well known that linear theory is totally 
inadequate for the analysis of stead^^ two-dimensional flows 
about airfoils when the Mach number approaches unity, or, 
in the present notation, when Jo approaches zero. 




Figure 34. — Variation of Cp with for various xjc. 

Ill oi-der to provide further information regarding the sig- 
nificance of the indicated variations of Cp with Jo, lines of 
constant J have also been included in figure 34. The local 
velocities are subsonic if J is negative and are supersonic if J 
is positive. The Mach number freeze is indicated in figure 
34 by the curve representing the variation of (\ with Jo at a 
given station xjc becoming parallel to a line of constant J. 
It can be seen from this figure that such behavior occurs at all 
stations on the airfoil for Jo near zero. 

Pressure drag. — The foregoing paragraplis have been 
concerned with the determination of pressure distributions 
on thin nonlifting circular-arc airfoil sections. Once these 
results are known, it is a simple matter to determine the 
section pressure drag coefficient Ca 



(c) Shock wave at trailing edge. 
Figuri: 33. — Concluded. 


( 103 ) 
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It is convenient at tliis point to introduce a I’edneed section 
drag coefTicient q defined in terms of reduced cpiantilies, 
thus 


Since tJie Mach man her freeze corresponds, in tlie approxi- 
mate theory to eciiiation (106) may be 

simj)lified to tlie following form: 


Ca 



ax 


(104) 


It is clear from the definitions of r, and x that tlie 
relation between Ca and Cd is 


(tic) 



^ dx 


(105) 


It should be remarked that some additional error is incurred 
in the present calculations of drag because tlie pressures are 
large and only i)oorly determined in the vicinity of the lead- 
ing and trailing edges. 

The variation of with has been computed and the 
results are shown in figure 35. It can be seen that is 
zero for less than —1.418 (corresponding to the critical 
Mach number). The rapid rise of q as is increased be- 
yond the ciitical is associated with the rearward movement 
of the shock and terminates abruptly when the shock reaches 
the trailing edge at d.825. The drag coedicient con- 

tinues to increase slowly with further increases of although 
at a much reduced rate; aiid, fiiially, at a value of some- 
what less than zero, ch becomes invariant with further changes 
in ^ 0 * This latter behavior is associated with the Mach 
number freeze, thus. 



(106) 



where (Zu)lb and (Zu)te refer to the ordinates of the upper 
surface at the leading edge and trailing edge, respectivel 3 ^ 
Since both of these quantities are zero for the circular-arc 
airfoils treated herein, it follows that 



The corres])onding results for symmetrical double-wedge 
profiles given in references 17, 13, and 14 by Trilling, Guder- 
ley and Yoshihara, and Vincenti and Wagoner, respectively, 
are also shown in figure 35. Exce])t for the value of ^ at 
^ 0 = 0 , tJie two sets of residts bear only a qualitative resem- 
blance. In some cases, the reason for the difference is 
clear; in other cases, the reasons are more obscure. An 
exami)le of the former concerns the residt that the drag of 
the circular-arc section is zero for less than —1.408, 
whereas the drag of the double-wedge section remains finite 
for all Tbe latter behavior results from the fact that the 
critical Mach number is zero for tlie double-wedge section. 
The reasons for tlie pronounced difference in shape between 
the two curves for drag ai subsonic speeds are not so clear. 
The immediate explanation is that the shock wave moves 
rearward across the chord of the circular-arc airfoil at a more 
rapid rate than it does with the double-wedge airfoil. For 
the circular-arc airfoil, the shock wave has moved to the 



Figure 35. — Theoretical drag results for circular arc and double wedge sections. 
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trailing edge when equals —0.825. This condition marks 
the end of the rapid increase in drag. On the other hand, 
Trilling’s results for the double-wedge airfoil indicate that 
the shock wave does not reach the trailing edge until the 
free-stream Mach number is unity, or Whether this 

difference is an actual proj)erty of the solutions of transonic 
small-disturbance theory for these two profiles, or the result 
of simplifying assumptions introduced in either the present 
analysis or that of Trilling remains an unanswered question 
at the present time. 

RESULTS IN TERMS OF PHYSICAL QUANTITIES 

The preceding section has summarized the results of the 
present calculations of the pressure distribution and drag of 
thin circular-arc airfoils in flows having free-stream Mach 
numbers up to unity. These results are given in terms of 
the reduced parameters Cp, c^, and which possess the ad- 
vantage of condensing the information for all thickness ratios 
onto a single curve, but the disadvantage of being some- 
what complicated and unfamiliar. Consequently, it is the 
aim of this section to re-express these results in terms of the 
more conventional quantities Cp, Ca, Mo, and t/c. 

Critical Mach number. — An airfoil property that is always 
of interest is the critical Mach number, Mcr- The variation 
of the critical Mach number with thickness ratio can be 
readily determined from the result that ^cr= — 1-408, thus. 

In this and in the remainder of the discussion, it is assumed 
whenever the residts of numerical computations are pre- 
sented that k is as defined in equation (7) and that y is 
equal to 1.4. Figure 36 shows a plot of the results of these 
calculations together with the corresponding results ob- 



Fkujre 3G. — Variation of critical Mach number with thickness ratio 
for circular-arc airfoil. 


tained using (a) linear theory and (b) thin-aiidoil theory for 
incompressible flow together with the Karman-Tsien rule to 
accoTint for the effect of compressibility. In both of the 
latter calculations, the variation of Cp^^ with Mo was deter- 
mined using equation (23), as is customary in engineering 
practice. It can be seen that the present calculations indi- 
cate that the critical Mach number is somewhat less than 
given hy either linear theory or by the Karman-Tsien method. 
This is in accord with the results found when the more exact 
theories are applied to thin airfoils. 

Pressure distribution. — The foregoing general results have 
been applied to determine the pressure distributions on a 
4-percent-thick circular-arc profile. These results are 
presented in two alternative forms in figures 37 and 38. The 
first of these shows the chordwise variation of Cp for various 
free-stream Mach numbers and illustrates the development 
and rearward movement of the shock wave as the free-stream 
Mach number increases beyond the critical. The second 
form shows the variation of Cp with Mo for selected points 
x/c on the airfoil chord. Three sets of auxiliary lines are 
also shown; they are (a) lines of constant local Mach number 
computed using equations (24) and (25) and (b) lines showing 
the variation of Cp with Mo at selected points on the airfoil 
chord computed using thin-airfoil theory for incompressible 
flow together with the Karman-Tsien rule. At Mach num- 
bers less than the critical, it can be seen that the pressure 
coefficient at a given point on the airfoil surface varies with 
Mach number in a manner similar to that predicted by the 
Karman-Tsien rule. It is also apparent, as was pointed out 
earlier in connection with equations (22) and (23) and figure 
2, that the value to be taken for the critical Mach number 
differs slightly depending on whether the lines of constant 
local Mach number are comi)uted by the present approxi- 
mation (eq. (24)) or by the exact isen tropic relation (eq. 
(25)). The irregular behavior of the pressures on the rear 
half of the airfoil at supercritical Mach numbers is associated 
with the passage of the shock wave. At Mach numbers near 
unity, a completely difl‘erent behavior is evidenced in which 
the values of Cp change with in such a manner that the 
local Mach number remains essentially constant. Direct 
calculation shows that the present theory implies the fol- 
lowing approximate relation for d(\l(lMo at a free-stream 
Mach number of unity.^ 



2 


/.*=.w(^ti) (110) 


This may be compared with the exact relation given by 
Vincenti and Wagoner in reference 14. 

G3^).u„=i^7+T“7+T 

* It might bo noted here that numerous investigators have given the appro.ximate relation 
/^\ 

y+l 

rather than that given in equation (110). As indicated by the auxiliary equations, the rela- 
tions for dCptdMo at Af»= 1 follow directly from the particular definition of k. It is interesting 
to observe that the expression for k given in equation (7) provides the more accurate relation 
for (dCp/dA/<.)A/o-i, even though, as noted in connection with the discussion of equation (96), 
it provides the lesser accuracy in the determination of (dMldMo)\fo~i- 
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Figure 37. — Chordwi.se variation of Cp for varoiis free-streain Mach numbers, 4-percent-thick circular-arc airfoil. 


Pressure drag. — The variation of the section pressure drag 
coefficient Ca with Afo has also been calculated for the 4- 
percent-thick circular-arc airfoil and the results are shown in 
figure 39. The general features of this curve are very 
similar to those discussed previously in connection with the 
correspondii'g curve of figure 35 for the reduced quantities. 
The major point of difference concerns the slope of the curve 
at a Mach number of unity. Thus, the preceding discussion 
has disclosed that dcdjdio is zero at a Mach number of unity; 
whereas figure 39 shows that dcaldAIo is negative at AIo=\. 
Direct substitution gives the following value for this slope: 



Vincenti and AVagoner have shown in reference 14 that the 
exact relation for flow about a closed airfoil is 



(113) 


The negative value of the slope given in equation (112) 
arises from tlie fact that the quantity k which appears in all 
the reduced parameters (q, (\, etc.) is a function of AIq. 
If, as in many other papers, k is equated to jUo and 

is thus independent of 71 /p, the value of the slope is zero, 
thus 



(114) 


COMPARISON WITH EXPERIMENT 

Inasmuch as the results (h'scribed in the preceding sections 
were calculated after making numerous simjilifications and 
approximations, not the least of which is the assumption of 
inviscid flow, it is desirable to include some comparisons 
with experiment. This is jiarticularly true for the present 
problem since it is well known that phenomena outside the 
scope of potential theory, such as separation, lioundary-layer 
shock-wave interaction, etc. are prominent features of tran- 
sonic flow about airfoils. Since it is indicated in the pre- 
ceding discussion that the present results are in general 
accord with the proven Prandtl-Glauert and Karman-Tsien 
results in the subcritical Mach number range, the following 
remarks will be confined to the supercritical range. 

Tliere are at least three papers available which present 
results of detailed measurements of flow at high subsonic 
velocities about symmetrical circular-arc airfoil sections, 
namely, reference 33 by Liepmann, reference 34 by Liep- 
mann, Ashkenas, and Cole, and reference 35 by AVood and 
Gooderum. The first two of these are concerned primarily 
with boundary-layer sliock-wavc interaction and contain 
statements casting doubt on the accuracy of the values in- 
dicated for the free-stream Mach number, a quantity of only 
secondary interest in their investigations. In general, it 
apiiears that values given for AIq in these two papers are 
somewhat too great. Tlu^ more recent investigation of 
AA'ood and Gooderum appears to be better in this particular 
and possesses the advantage of being made with an inter- 
ferometer so that knowledge is gained about the entire flow 
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Figure 38. — Variation of Cp with Mo for various values of x/c, 4-percent 
thick circular-arc airfoil. 

field. A disadvantage, however, is that the tests were made 
using a 12-percent-thick model, which severely strains the 
assumptions of transonic small-distin-bance tlieory. Never- 
theless, the comparisons will be made with the data of 
reference 35. 

The studies of Liepmann (ref. 33), Ackeret, Feldmann, and 
Rott (ref. 36), and others have shown that the boundary 
layer can have a profound influence on a transonic flow field. 
This immediately^ raises a question regarding the usefulness 
of a potential-flow theory, such as the present, which dis- 
regards the boundary layer completely. In order to illustrate 
better the nature of these effects, two interferograms of the 
flow about 12-percent-thick circular-arc airfoils are repro- 
duced from reference 35 and shown in figure 40. The free- 
stream Mach number is 0.88, and the Reynolds number based 
on the chord is 600,000 for both photographs. The condi- 
tions for the two flows differ in that the boundary layer is 
laminar in the flow pictured on the left and turbulent in 
that pictured on the riglit. The interferogram for the 
laminar case shows that the shock waves are of the X type 
and That the flow separates near the midchord station. 
When the boimdaiy layer is turbulent, however, it may be 
seen that the shock wave is of the simple single wave type, 



Mo 

Figure 39. — Variation of Cd with Mo, 4-percent-thick circular-arc 

airfoil. 

and that the extent of the region of separated flow is greatly^ 
diminished. From this pair of photographs, it is apparent 
that two of the simplifying assumptions introduced in the 
course of the present analy^sis (i. e., (a) no flow separation 
occurs and (b) the shock wave is a single normal wave) are 
in better accord with the physical phenomenon if the bound- 
ary lay-er is turbulent ahead of the shock wave than if it is 
laminar. Figure 41 shows experimental pressure distribu- 
tions determined from the two interferograms of figure 40 
together with the corresponding theoretical results. Because 
of difficulties in interpreting the interferograms, the experi- 
mental pressure distributions presented in reference 35 and 
reproduced in figure 41 are terminated at the separation 
point. It can be seen that the theoretical results are in 
substantial agreement with the experimental data available 
for the portion of the airfoil forward of the separation point. 
Although the experimental and theoretical pressure distri- 
butions matched equally^ well for all the test data of reference 
35, the degree of correspondence must be attributed, in part, 
to a fortunate cancellation of errors since individual com- 
ponents involved in the perturbation analvsis (e. g., the 
relation between w and dZjdx, Cp and the perturbation 
velocity components, etc.) contain appreciable errors when 
applied to such a thick airfoil. Although the corresponding 
experimental values are not availalile for the pressures at 
stations aft of the separation point, it is presumed from other 
experimental data that the pressures reach their maximum 
negative value in the vicinity of the separation point and 
return toward tiie free-stream value at stations farther aft 
on the airfoil, rather than varying in the manner indicated 
b.v the calculated results. 

Before leaving this topic, the following remarks should 
be made concerning experimental factors which may affect 
these conclusions. The first is that the authors of reference 
35 do not consider the data for the turbulent boundary- 
lay-er to be as reliable as that for the laminar boundary 
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(a) Laminar boundary layer. 

(b) Turbulent boundary layer. 

Fkh re 40. — Interferograins of flow about 12-pereent-thick circular-arc 
airfoil at a Mach number of 0.88 and a Reynolds number of 600,000. 
From reference 85. 

layer. The reason for this is that the method used to make 
the l)Oiindarv layer tuiTulent (the upper half of the airfoil 
was mounted on a flat ])late which extended one chord 
length forward of tlie leading edge of the airfoil) produced 
such a thick houndary layer that difficulties were encountered 
in correctly extrapolating the lines of the interferograms to 


tlie airfoil surface.'^ The differences in boundary-layer 
thickness are clearly evident in the interferograms shown 
in figure 40. The second stems from the fact that the tests 
were conducted with a 12-i)ercent-thick airfoil. Compari- 
sons of theory and ex])eriment for such thick airfoils not 
only strain the small-disturbance assumptions of the theoiy 
but also em])hasize unduh^ such features of the flow as the 
curvature of the shock wave and ])ound ary-layer se]3aration 
which are disregarded coni])letely in the theory. The third 
is concerned with the fact that results of a recent flight 
investigation by Ilarrin (ref. 37) have shown that, at Rey- 
nolds numbers of the order of 20,000,000, there is very 
little difference in ])ressure distribution with laminar and 
turbulent boundary layers. In particular, it is found that 
the X type shock and large region of separated flow commonly 
associated with laminar boundary layers do not occur. If 
this result is substantiated by further experiments, it will 
l)e of particular significance in work such as the ])resent 
where some assum])tion has to be made about the nature 
of the shock system, and the assumption that the shock 
wave is a single wave leads to the greatest simplification. 

The principal discrepancy l)ctween calculated pressui’e 
distrif)Utions and those measured with a turbulent ])oundary 
layer stems from the fact tliat the shock wave meets the 
airfoil surface at a station farther forward than calculated. 
44ie interferograms of figure 40 show that a substantial 
portion of the forward shift is due to a pronounced curvature 
of the shock wave near the airfoil. This suggests a comj)ari- 
son of the calculated shock ])ositions with not only the 
positions observed at the airfoil surface but also at some 
distance, say a half chord length, away from the airfoil. 
The results of such a comparison with the data for turbulent 
l)Oundary layers given in reference 35 are shown in figure 42. 
It can be seen that the calculated shock positions vary with 
Mo in a manner which i)arallels the experimental results at 
= but not those at the airfoil surface. It is almost 
certain that the details of this ])henomenon are greatly 
influenced l)y separation and by interaction of the shock 
wave and boundary layer, and are, hence, beyond the reach 
of separation-free ])otential theory for transonic flow. 
Since the pressure gradients and shock strength increase with 
airfoil thickness, it might be j)resunied that these effects would 
be greater for thick airfoils than for thin airfoils. An indi- 
cation of such a trend is furnished by the pressure-distribution 
data of Gothert (ref. 38) for XACA 0006, 0009, 0012, 0015, 
and 0018 airfoils. Figure 43 shows a plot of the variation 
with ^0 (^f flic shock position at the surface of each of these 
airfoils. If the shock positions were changing in accord 
with the similarity rules of transonic flow theory, these 
results would all fall on a single curve. It can be seen, 
however, that this is not the case and that the shock moves 
rearward across the chord at a slower rate for the thick 
aii’foils than for the thin ones. It appears, therefore, that at 
best, the calculated shock positions will only agree with those 
found experimentally for very thin airfoils. 

® Results are also given in reference 35, although not for a Mach number of 0.88, in which 
the boundary layer is made turbulent by a wire stretched across the test section one chord 
length ahead of the leading edge. When these results are plotted in the form shown in figure 
39, they fall about half way between the calculated results and those obtained with the half 
model mounted on a plate. 


Shock position, x/c 
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(a) Laminar boundary layer. (b) Turbulent boundary layer. 

Figure 41.— Comparison of theoretical and experimental chordwise pressure distributions for a I2-percent-thick circular-arc airfoil at a :\Iach 

number of 0.88. Experimental data from reference 35 for a Reynolds number of 600,000. 


Theory Experiment 

(Turbulent flow) 



Figure 42.— Variation of shock position with Mach number for 12- 
percent-thick circular-arc airfoil. 



Figure 43. — Effect of thickness on shock position. 
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It is apparent that ciiaiiges in the shock position will he 
accompanied, in many cases, by substantial changes in the 
])rcssure drag of the airfoil. In particular, the forward shift 
of the shock noted for the thick airfoils will diminish the 
region of high negative pressures on the rear of the airfoil, 
thereby decreasing the drag. These effects will, of course, 
be zero until the critical Mach number is exceeded, will 
increase as the shock moves across the airfoil, and may 
diminish somewhat again as sonic free-stream velocity is 
approached, since the actual, as well as calculated, shock 
position then approaches the trailing edge. These effects 
would probably alter the curve for the pressure drag of a 
thick airfoil so that it would appear somewliat like that 
illustrated in figure 44. 



Figure 44. — Presumod drag results for a thick airfoil. 


CONCLUDING REMARKS 

The foregoing results are encouraging in that they show 
that the introduction of a small number of rough, although 
reasonable, assumptions leads to a relatively simple method 



(a) Theory. 

Figure 45. — Comparison of theoretical and experimental lines of constii 
0.88. Experimental data from reference 35 for Ueyi 


for the calculation of pressure distriliutions on thin circular- 
arc airfoils at all Mach numbers up to unity. Perhaps the 
most important aspect of the present work is the recognition 
of the quadratic imture of the integral equation and the 
retention of this feature in the iteration solution. With the 
knowledge that acceptable results can be obtained without 
excessive effort, it appears worthvvdiile to re-examine the 
approximate solution of the equations with an eye toward 
im])rovement, or elimination, of the simplifying assumptions. 
Proba])ly one of tiie weakest elements of the present method 
is the velocity-distribution function introduced in equation 
(63). Tliis particular function is used to determine the 
entire flow field but, actually, only insures that the velocity 
and velocity gradient have the correct value at the airfoil 
surface and that the velocity diminishes toward zero at 
infinity as Ifz^. As can he seen by comparing the experi- 
mental and calculated results shown in figure 45, this func- 
tion succeeds reasonably well for circular-arc airfoils. On 
the other liand, such a simple function cannot be expected 
to give good results for all airfoils. For instance, the present 
velocity-distribution function cannot be expected to provide 
good results for airfoils having flat surfaces over a substantial 
part of the chord (e. g., wedge airfoils, etc.) since it indicates 
no attenuation with distance above an element of the airfoil 
where the radius of curvature of the surface is infinite. Gull- 
strand (ref. 22) has proposed a different velocity-distribution 
function which satisfies one more known requirement, but 
it lias not been established as yet whether or not it is suffi- 
ciently general to cover all interesting cases. 

Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Nov. 19 y 1953 



(b) Experiment. 

Cp around a 1 2-percent-thick circular-arc airfoil at a Mach number of 
Is number of 600,000. Turbulent boundary layer. 
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PRINCIPAL SYMBOLS 


speed of sound 

speed of sound in tlie free stream 

critical speed of sound 

function defined in equations (63) and (64) 

pressure coefficient, ^ — — 

Pott 2 

cy ^ 0 

(Uoky^^ ^ 

(f/cY'^ ” 

chord 

section pressure drag coefficient, 




■Cd 


(Uokr 

pressure drag 

function defined in equation (67) 
function defined in equation (71) 
function defined in equations (72) and (73) 
function defined in equations (76) and (77) 
function defined in equations (57) and (58) 
coefficient of nonlinear term of differential equation 
for ip (wSee eqs. (7), (17), (18), and (21).) 

2 Ul—\ 

width of element used in approximating the 
chordvvise velocity distribution 
local Mach number 
critical Mach number 
free-stream Mach number 
inward normal to surface enclosing Ti 
arbitrary point cr, 2/, 2 
static pressure 
free-stream static pressure 
region of integration 

^ {x-lY+(z-fy 

^^(x-^y+{y~vy+i^-kr 

surface of shock wave 
wing semispan 

maximum thickness of profile 
free-stream velocity 

perturbation velocity components parallel to 
X, y, z axes, respective!}^ 
k 

k 


k 


u, ?, w 

Cartesian velocit}^ components on the two sides of 
a shock wave, with u being parallel to the flow 
direction ahead of the shock 


I— f 


b' 

2^; 2/, ^ 

Cartesian coordinates where x extends in the di- 
rection of the free-stream velocity 

X 

X 

y 


z 


Z 

ordinates of wing profiles 

Z 

kilo y 

0 

V17-W 

7 

ratio of specific heats, for air 7=1.4 

A 

difference between values of quantity on the upper 
and lower sides of the xy plane 

^0 

variables of integration corresponding to x^y^z 
l-Mo^ 

l-l\p 

[Uokmr^^ 

Po 

free-stream density of air 

s 

surface enclosing volume R 


kUo t 

T 

c 

<!> 

velocity potential 

'P 

perturbation velocity potential 


k 



'p 

1 /W\ 

^’3 dp \ 2 / 

Subscripts 

a 

values ahead of shock wave 

b 

values behind shock wave 

cr 

conditions associated with the critical Mach 
number 

L 

values given by linear theory 

LE 

leading edge 

1 

values on the lower surface of wing or wake 

A/„=l 

values at 

P 

values at arbitrary point P 

S 

values along shock wave 

TE 

trailing edge 

u 

values on the upper surface of wing or wake 

W 

values at the wing surface 
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DISCUSSION OF CERTAIN PROPERTIES OF THE QUADRATIC ITERATION METHOD 


INTRODUCTION 

The present results are obtained using an iteration process 
wliich differs from that of previous workers in that the 
quadratic nature of the integral eciuation is recognized and 
retained throughout the analysis. This apiiendix is con- 
cerned with an investigation of certain properties of the 
two methods and a discussion of the applicability of eacli 
method to transonic-flow problems. Certain other problems 
will be avoided by confining attention to symmetrical non- 
lifting wings with sufficiently smooth entries at leading and 
trailing edges that there arc no stagnation points. Methods 
of handling stagnation points are discussed in references 39 
and 40. 

Tlie fundamental relations for the following discussion, 
equations (41) and (43), can be written in the following more 
abbreviated form 

U = Tir.-h (Bl) 

»=S.+|-t (B2) 

where 7i and h refer to the integrals o\cv U and S in equa- 
tions (41) and (43), respectively, and Ul is the linear-theory 
solution given by equation (42). As is evident from the 
derivation in the text, equations (Bl) and (B2) are totally 
equivalent since the latter is derived from the former by 
partial int(‘gration such that 

A=f-| m 

Since metJiods for inverting, or solving, nordinear integral 
equations such as equations (Bl) and (B2) are not known, 
tlie present analysis, in common witli many pz-evious analy- 
ses, seeks an apizroximate solution by iteration. These 
methods pi’oceed, in general, by substituting a knowji 
function for u into the integrals, integrating, and solving the 
resulting algel)i*aic equation for a new appi*oximation for z/, 
which is in turn substituted into the integrals, etc. The 
pi’Ocedui*e eznployed in the pi*esent analysis differs from that 
traditionally employed in that it is based on equation (B2) 
rather than equation (Bl). The most obvious difference 
between tlie two procedures is that the algebraic equation to 
solve for the new u is linear in the traditional analysis and 
quadratic in the jiresent analysis. Although solutions of 
equations (Bl) and (B2) must actually be identical, it is 
sliown in equations (91) through (94) that this is not always 
the case for the results indicated by a partially completed 
iteration calculation. Since the linear iteration method is 
accepted universally for purely subsonic flow, and by some 
for mixed flows, it is im])ortant that we examine the relation 
between (he results of the two procedures in greater detail. 


SHOCK-FREE FLOW 

Shock-free flow will l)e considered first. The integrals 7i 
and I 2 are then simply 

R 

'““-[sJJJ I i(i;) 


and will l)e designated 7ifiZ-) and to call attention to 

the fact that they involve tlie square of zZ. Tlie linear-llieory 
solution Ul has the form 

Ul=~tQi(x, y, i) = ui7 (B6) 

where ai is a known function. Both the linear and quadratic 
iteration procedures can be started by letting u = 0 and desig- 
nating this approximation by uq. 

The results for successive iteration steps using the linear 
iteration method are determined by solving 

Un-\-l^=Ur~ I\(Un) (B7) 

In this way we obtain 


ll\^UL = (hT 

U 2 = Ul— I\ (u\^) == Til — r-7, (a 1 2) = ai 7 + C 2 
u-i = Vl — I\ (z/2^) = c , 7 + + 0 (? ^ 


(B8) 

n = l 

where an=an(Xj y, 2 ) are functions of the wing geometry, 
and ultimately therefore of Oi. This procedure has the 
interesting property that one and only one additional term 
in the power series expansion for u is determined in its final 
form with each iteration step. That is 

^/n-i + ^ «T" + 0(r''''^^) (B9) 

Note that a single expression is obtained for all Mach num- 
bers and thickness ratios, and that the magnitudes of these 
parameters are reflected solely in the value for r. The 
principal question concerns the range of convergence of the 
series. It is evident that the series cannot continue to con- 
verge as Mo approaches unity because r approaches infinity, 
but tlie precise liound on the radius of convergence is difficult 
to establish at this point. 
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The results for successive iteration steps using the quad- 
ratic method are determined by solving 


Thus 

where 

and 


^n+1 9 


Un+1—^ ± A^i 


— l±'\fl — Xn 

Xn = 2UL — l2{Ur?) 

x„<l 


(BIO) 


(Bll) 


if Un^\ is to remain real. The minus sign is associated with 
local velocities that are subsonic and the plus sign with those 
that are supersonic. This result can be made more nearly 
comparable with the power series representation provided 
by the linear iteration method b}^ expanding the radical of 
equation (Bll) under the assumption that |xn|<Cl> whereby 
we obtain 


t/n+l = l± 




X 


2 , 

n 




Xn^— 



Consider first purely subsonic flow. Then successive 
iterations 3 ueld tlie following results 


FLOW WITH SHOCK WAVES 

Although many of the general notions described for shock- 
free flow carry over into the analysis of transonic flows with 
shock waves, certain changes must be introduced. It is 
clear, for instance, that the linear iteration method cannot 
be used with the starting u distribution taken as zero or ul, 
because then equation (B8) would be reproduced and no 
discontinuities or sliock waves would ever appear. On the 
other hand, it miglit seem reasonable to suppose that a 
successful iteration calculation could be accomplished using 
a starting u distribution containing a shock wave and esti- 
mated to be close to the correct solution. This is, in fact, 
what is done in the present analysis, but using the quadratic, 
rather than linear, iteration method. Although, as showm 
in the preceding section, there is no essential difference be- 
tween the results of the two iteration methods in the sub- 
critical range, the quadratic method emerges as definitely 
superior when the free-stream Mach number approaches 
unity. Its decisive advantage lies in the fact that small 
errors in the estimated u distribution lead only to finite errors 
in the calculated u as approaches unit}^ whereas infinite 
errors can result from the calculations performed by the 
linear iteration method. The means by which small errors 
can be magnified in the linear iteration process can be seen 
by examining equation (Bl) together with the definition of 
u and ul- 

Jc Jc 

u; 


« 0=0 


Tims equation (Bl) can be rewTitten as follows; 


Ui=UL+0(UL^=air+0(T^ 


= CilT 


(B13) 

n=l 

wliere the an have the same values as in equation (B8). 
This series and that given by equation (B8) are the same 
since they can be made to agree to any arbitrary number of 
terms. Thus the two methods lead to the same result for 
purely subsonic flow. AVhereas there is some doubt about 
the range of convergence of equation (B8), it is clear that 
equation (Bl3) only applies to purely subsonic flow. This 
is evident, first of all, from the use of the minus sign in equa- 
tion (Bll). Moreover, since the series only converges for 
lxnl<Clj 111® largest value for u which can be represented by 
equation (Bl3) is unity. Since u equals unity when the local 
velocity is sonic, this result can be interpreted as indicating 
that equation (B8) converges for purel}^ subsonic flow, but 
diverges for mixed flow. 

There also exists a range in which the expression given in 
equation (Bl2) applies using the plus sign for certain regions 
of the flow field. The resulting series expansion for u, how- 
ever, does not lead to that given by equation (B8) or (Bl3), 
inasmuch as the first term is independent of 7. 


u = u^—ll 



Now, as Mo approaches unity, u remains finite but Ul ap- 
proaches infinity. The largest values of are attained on 
the airfoil surface where approaches infinity as 1//3 in 
two-dimensional flow and In in three-dimensional flow. 
This means, in the linear iteration method, that the desired 
solution u must be determined from the small difference be- 
tween two large numbers. Since 7i is only evaluated ap- 
proximately in an iteration method by substituting a function 
for u which differs from the exact function by a finite amount, 
it is apparent that very large or infinite errors may result as 
Mo approaches unity. On the other hand, equation (B2), 
which is used in the quadratic iteration method, can be 
rewritten as follows: 


where 


0 h 



Now in this case the term containing Ul becomes indefinitely 
small compared with as the Mach number approaches 
unity, and only finite errors result when an approximation 
is substituted for u in lo. 
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